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Abstract 

In Part I of this series of papers we have described a general formalism 
to compute the vacuum effects of a scalar field via local (or global) zeta reg¬ 
ularization. In the present Part II we exemplify the general formalism in a 
number of cases which can be solved explicitly by analytical means. More in 
detail we deal with configurations involving parallel or perpendicular planes 
and we also discuss the case of a three-dimensional wedge. 
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1 Introduction 


This is the second part of our series of papers about zeta regularization and vacuum 
effects for a scalar field. In Part I [19] we have considered a neutral, scalar quantum 
field on (d + l)-dimensional Minkowski spacetime, assuming the field to be confined 
within a spatial domain and to fulfill suitable boundary conditions; we have also 
indicated the possibility to replace hi with a Riemannian manifold, or either an open 
subset of it (and the Minkowski environment with a curved ultrastatic spacetime). In 
the same work we have been presenting general methods for the zeta regularization 
of the vacuum expectation value (VEV) of the stress-energy tensor, of the total 
energy and of the boundary forces, emphasizing the fact that we were producing 
a set of general rules to be applied almost mechanically in specific configurations. 
The illustration of such mechanical rules was started in Part I with some simple 
examples in one spatial dimension, and it is continued in the present Part II with 
more engaging configurations. Here we consider a number of cases, in which the 
necessary computations can be carried over by purely analytical means; this makes a 
difference with respect to the subsequent Parts III [20] and IV [21], where application 
of the general rules will require a mixture of analytical and numerical methods. 
Most of the completely solvable cases considered here have been previously treated 
in the literature, typically in the case of spatial dimension three, often with special 
choices of the boundary conditions and, in most cases, considering only the conformal 
part of the stress-energy VEV; for each one of the cases already treated a specific 
approach has been employed, possibly different from zeta regularization. 

The present work attains more generality as for the spatial dimension, the bound¬ 
ary conditions and the presence of a non-conformal part in the stress-energy VEV; 
the second, perhaps more significant, contribution of this paper is the unified view¬ 
point mentioned in the previous comments, that is, the presentation of all cases as 
applications of the same apparatus. 

To be more specific, the configurations that we analyse are the following, 
i) First of all we consider a massless field in odd spatial dimension confined between 
two parallel (hyper-)planes, for several kinds of boundary conditions; see Section 
[3] We have already considered this model in spatial dimension d = 3 for Dirich- 
let boundary conditions in HU; therein, the analytic continuation required by zeta 
regularization was performed using ad hoc, known results on the special functions 
related to this specific configuration (namely, the Riemann zeta and the polyloga- 
rithm). As a matter of fact, the literature on the configuration with two parallel 
planes is immense, both regarding local and global aspects; here we only cite a few 
references. In his seminal paper [9], using an exponential cut-off regularization along 
with Abel-Plana resummation, Casimir was the first to compute the total energy 
and the boundary forces for the case of two parallel planes; concerning local aspects, 
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the foremost derivation of the full stress-energy tensor VEV was given by Brown and 
Maclay jSj, using a point-splitting technique (jl). Computation of both global and 
local quantities for this model was later reproposed by several authors, using various 
regularization techniques: see, e.g., the monographies by Milton [25], Elizalde et ah 
PUDS], Bordag et ah [3] (see, as well, the works cited therein) and the papers by 
Zimerman et ah [3D] and by Esposito et ah [TB]. We will give further bibliographical 
references in Section [3] In the present paper we resort to the formalism of integral 
kernels developed in Part I in order to derive automatically the analytic continua¬ 
tions required by the use of zeta regularization. 

ii) Next, we consider a massive held between an arbitrary number of perpendicular 
(hyper-)planes, fulfilling either Dirichlet or Neumann boundary conditions on each 
one of them; see Section [I] To the best of our knowledge, this type of configuration 
was only considered by Actor [I] and by Actor and Bender [2]; both these papers 
deal with a scalar held fulfilling Dirichlet boundary conditions on at most three per¬ 
pendicular planes (along with other similar models in three spatial dimensions, all 
involving boundaries consisting of hat, perpendicular parts). More precisely, in [1] 
the author hxes d — 3 and evaluates the renormalized effective Lagrangian, along 
with the VEV (O|0 2 (x)|O) ren (plus their analogues for non-zero temperature); in [2] 
d is arbitrary and it is computed the VEV of each component of the stress-energy 
tensor. In both works cited above, a zeta type approach based on the use of heat ker¬ 
nels is employed; however, contrary to our regularization scheme, this approach also 
involves the subtraction of terms (corresponding essentially to Minkowski spacetime 
contributions) which diverge for any value of the regulator parameter. Besides, let 
us stress once more that our method is more systematic and the construction of the 
required analytic continuations descends automatically from a general framework. 

iii) Finally, we consider a massless held confined within a wedge of arbitrary width 
in spatial dimension d = 3, for several types of boundary conditions; see Section 
O We also consider a variation of this configuration, which corresponds essentially 
to identify the sides of the wedge; this is the so-called case of the “cosmic string” 
(see subsection 15.91) . Some of these cases have already been treated by Dowker et 
ah [12, 13], Deutsch and Candelas [Tl| (also discussing the electromagnetic case) 
and, more recently, by Saharian et ah [29], 31] and by Fulling et al. [S3] (see also 
the citations in these works and [SHUM ; nearly all of these authors use the point 
splitting approach, or some variant of it. More in detail, in jT2j and m attention 
is restricted to the conformal part of the stress-energy VEV for either Dirichlet or 
Neumann boundary conditions, while in [29] [3T] also the non-conformal part is con¬ 
sidered, but in the Dirichlet case only; in [23], instead, the authors only show the 
graphs of the energy density and of the pressure components (for which no explicit 

2 Actually, both Casimir and Brown-Maclay considered the case of electromagnetic field; yet the 
methods employed by these authors can be trivially adapted to the case of a scalar field. 
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expression is given), derived via a point-splitting approach for several configurations 
and various choices of the parameters describing the theory. Our approach via zeta 
regularization considers several types of boundary conditions, both in the conformal 
and in the non-conformal case; in the subcases already analysed in the literature 
cited above, we obtain the same results. 

Before describing the applications (i-iii), in the forthcoming Section [2] we summarize 
the general scheme of Part I (more precisely, the parts of it required by the appli¬ 
cations analysed in the present paper). This summary has been written just for the 
comfort of the reader who, in absence of it, would be forced to skip continuously to 
Part I to recover the basic identities on integral kernels and analytic continuation 
employed here. 

Finally, we point out that some of the computations presented in this paper have 
been performed using Matheraatica in the symbolic mode. 

2 A summary of results from Part I 

2.1 General setting. Throughout the paper we use natural units, so that 

c — 1 , h= 1. (2.1) 

Our approach works in (d+ l)-dimensional Minkowski spacetime, which is identified 
with R d+1 using a set of inertial coordinates 

x = (x M ) /1= o,i....,«i = Or°,x) = (t,x) ; (2.2) 

the Minkowski metric is (r ] fiu ) = diag(—1,1,..., 1). We fix a spatial domain 0 C R c/ 
and a background static potential V : O —> R. We consider a quantized neutral, scalar 
field (j) : R xO —» £ sa (i?) (5 is the Fock space and £ sa ($) are the selfadjoint operators 
on it); suitable boundary conditions are prescribed on <90. The field equation reads 

0 = {-d u + A - V(x))?(x,t) (2.3) 

(A := Y2i=i is th e d-dimensional Laplacian). We put 

A:=-A + V, (2.4) 

keeping into account the boundary conditions on <90, and consider the Hilbert space 
L 2 (0) with inner product (f\g) f Q dx f(x)g(x). We assume A to be selfadjoint 
in L 2 (0) and strictly positive (i.e., with spectrum cr(A) C [e 2 , +oo) for some e > 0); 
the latter assumption is sometimes relaxed requiring only that A is non-negative 
{cr{A) C [0, Too)). 
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We often refer to a complete orthonormal set (F k ) k ^ of (proper or improper) eigen¬ 
functions of A with eigenvalues (uiDkeic (<^>k ^ £ for all k G /C; with the relaxed 
condition of non-negativity, we only have uj k ^ 0). Thus 

F k : D —> C; AF k = u k F k ; 

= 7(7, h) for all k, h G K. . 

The labels k G /C can include both discrete and continuous parameters; dk indi¬ 
cates summation over all labels and 5(h, k) is the Dirac delta function on fC. 

We expand the field 0 in terms of destruction and creation operators correspond¬ 
ing to the above eigenfunctions, and assume the canonical commutation relations; 
|0) G S' is the vacuum state and, as already indicated, VEV stands for “vacuum 
expectation value”. 

The quantized stress-energy tensor reads (£ G R is a parameter) 

%u ■■= (1 - 2£) d^cj) o d u (f> - Q - 2£^J riuJd^dxd) + V0 2 ) - 2£ 0 o 0; (2.6) 

in the above we put A o B := (1/2 ){AB + BA) for all A, B G £ sa (5"), and all the 
bilinear terms in the field are evaluated on the diagonal (e.g., C00 o d v (j) indicates 
the map x i—>■ d^ix) o d v <f>(x)). The VEV (O|T0,|O) is typically divergent. 

2.2 Zeta regularization. The zeta-regularized field operator is 

0“ := (k- 2 ^)- u/4 J , (2.7) 

where A is the operator (12.41) . u G C and k > 0 is a “mass scale” parameter; note 
that 4> u \ u= o = 0, at least formally. The zeta regularized stress-energy tensor is 

f; v := (l-20^? u o^?“-Q-2^//^(9 A M? u + n? u ) 2 ) -2£0“o«V0“. (2.8) 

The VEV (0|T“J0) is well defined for large enough (see the forthcoming subsec¬ 
tion [23]); moreover, in the region of definition it is an analytic function of u. The 
same can be said of many related observables (including global objects, such as the 
total energy VEV). 

For any one of these observables, let us denote with F{u) its zeta-regularized version 
and assume this to be analytic for u in a suitable domain WoCC. The zeta approach 
to renormalization can be formulated in two versions. 

i) Restricted version. Assume the map W 0 G C, u (->• F(u) to admit an analytic 
continuation (indicated with the same notation) to an open subset U C C with 
7/30; then we define the renormalized observables as 

Fren ■= F(0) . (2.9) 
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ii) Extended version. Assume that there exists an open subset U C C with Uq <ZU , 
such that OgW and the map uEUq (->• F{u) has an analytic continuation to 7/\{0} 
(still denoted with F). Starting from the Laurent expansion F{u) = Ylt=-oo J r k'u k , 
we introduce the regular part (RP J 7 ) (u) := Y^k=Q^kU k and define 

F ren := (RP F)(0) . (2.10) 


Of course, if T is regular at u = 0 the defnitions (12 .OR (12.1 Oft coincide. 

In the case of the stress-energy VEV, the prescriptions (i) and (ii) read, respectively, 


(0|f„„(i)|0)„„ := (0|f“(U|0) , 


( 2 . 11 ) 


(0\T^(x)\0) ren RP 


o (0|f" (x)|0) . 


( 2 . 12 ) 


2.3 Conformal and non-conformal parts of the stress-energy VEV. These 
are indicated by the superscripts (0) and (■), respectively; they are defined by 


(0|f^|0) ren = (0\T$>\0) ren + (£-&) (0|fW|0) ren , 


where we are considering for the parameter £ the critical value 


&:= 


d—1 

~AcT 


(2.13) 


(2.14) 


2.4 Integral kernels. If B is a linear operator in L 2 (fl), its integral kernel is the 
(generalized) function (x,y) 6 fixfl 4 £>(x,y) := (<5 x |£>(f y ) (5 X is the Dirac delta 
at x). The trace of B , assuming it exists, fulfills Tr£? = J Q dxB(x. x). 

In the following subsections A is a strictly positive selfadjoint operator in L 2 (fl), 
with a complete orthonormal set of eigenfunctions as in Eq. (12.5f) : in some situations 
(explicitly indicated) we only require A to be non-negative. In typical applications, 
A is the operator (12.4p . 

2.5 The Dirichlet kernel and its relations with the stress-energy VEV. 

For (suitable) s £ C, the s-th Dirichlet kernel of A is 

r f ir _ 

D s (x,y) := -4 _s (x,y) = / — F k {x)F k (y) . (2.15) 

Jk. u k 

If = —A + V (with V a smooth potential) is strictly positive, the map D s ( , ) : 

x —» C, (x, y) i —> D s (x, y) is continuous along with all its partial derivatives 
up to order j £ N, for all s £ C with JLs > d/2 + j/ 2; in particular, D s and its 
derivatives up to order j are continuous on the diagonal y = x, a fact of special 
interest for our purposes. (Under stronger assumptions on A, one can give results 
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of absolute and uniform convergence of the eigenfunction expansion in Eq. (I2.15I) . 
for suitable values of s; the same can be said for all the corresponding derivatives). 
Recalling Eq. (12.8p . the regularized stress-energy VEV can be expressed as follows: 


<OWo(x)|0) = «' 


1 (x, y) + (\-( ]Ur(J, ! , + Vix)tl), : i (x, y) 


, (2.16) 


y=x 


= K 


<0|T„"(x)|0) = (0|T“(x)|0) = 0. 
(0K“(x)|0) = (0|7™(x)|0) = 

- (3 I< S/ + ' / (x))£> a ±i(x,y)) + 

1 


(2.17) 


(2.18) 


dx i yj ^ 9 x i x j 1 Du±i (x, y) 


y=x 


((0|T“ ;/ (x)|0) is short for (0\T™ u (t, x)|0); indeed, the VEV does not depend on t). 

Of course, the map O —> C, x (0|T“ iy (x)|0) possesses the same regularity as 
the functions x G h H- Du±i (x, x), d zw Du±i (x, x) (z,w any two spatial variables); 

so, due to the previously mentioned results, x (->■ (0|T“ jy (x)|0) is continuous for 
3 R,u > d + 1. 

The renormalized stress-energy VEV is (0|T Atl/ (x)|0) ren := RP\ u=0 (0\T^(x.)\0); in¬ 
troducing the functions 


D±1 ( x ,y) := RP jK u D^(x,y)) , 


u =0 


d zw D[ K \-x,y) := RP (n u d zw Du± i(x,y)) 

2 u=0 V 2 / 

(with z,w any two spatial variables), this can be expressed as follows: 


(0|r 00 (x)|0) ren = 


i +p £><|(x, y) + (i-p (cf‘g v , + V (x))£|W(x, y) 


Mf 

2 


y— x 


(01 Tqj (x) 10) ren = (0|f iO (x)|0) ren = 0 , 
(0|^p( x )|0) ren = (0|T ii (x)|0) ren = 

- y) - (9 x ^+V(x))L>||(x,y)) + 


(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


+ ((2 - - £ d xi x i) D { + |(x, y) 


y=x 














If D «±i(x, y) and d zw Dt i+i(x, y) have analytic continuations regular at u — 0, indi¬ 
cated with _D ± i(x,y) and d zw D i(x, y), for any n > 0 one has 


D ±l ( x >y) = ^±i( x >y)» 

2 2 

dzwD i K) (x, y) = i (x, y) . 

2 2 


(2.24) 


In the sequel we will consider the total energy VEV and express it in terms of the 
trace Tr*4. s , fulfilling 

Tr A~ s — / dx.D s (x, x) . (2.25) 

Jn 

2.6 The heat and cylinder kernels. For t G [0, +oo), these are given by 


AT(t;x,y):=e U (x,y) = dk e tw * A fc (x)A fc (y) ; (2.26) 

Jk 

T(t;x,y) := e~ tva (x,y) = [ dk e~ tuJk F k (x)F k ( y) . (2.27) 

Jk. 

Sometimes we also consider the modified cylinder kernel 

T( t;x,y) := (VA V iV3 )(x,y) = f — e~ tuJk F k (x)F k ( y) ; (2.28) 

Jk u k 

it turns out that T(t; x, y) = — d{t (t; x, y). 

If A = —A + 1/(F smooth) is strictly positive, the map K( t; , ) : x —> C, 
(x, y) i—^ A'(t;x, y) is continuous along with all its partial derivatives of any order, 
for all t > 0 (the same holds for T and T ); with stronger assumptions on A, 
the eigenfunction expansions in Eq.s (I2.26jl (12.271) (12.281) converge absolutely and 
uniformly, along with all their derivatives. 

If A is not strictly positive, but non-negative, Eq.s (I2.26j) (I2.27j) continue to make 
sense; moreover, if 0 is in the continuous spectrum of A, its spectral measure vanishes 
and the operator \fA e~ tv ^ can still be defined; in this case, typically, the modified 
cylinder kernel and its eigenfunction expansion in (I2.28jl still make sense. 

The cylinder trace, if it exists, is 

T{ i)-=Tre- tVJ = [ rfxT(t;x,x). (2.29) 

Jn 
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2.7 The Dirichlet kernel as Mellin transform of the heat or cylinder ker¬ 
nel. For suitable values of sGC (see Part I), there hold 


^ r +°o 

D s(x,y) = ^TT dt P _1 K(i ; x, y) ; (2.30) 

r (s) Jo 

^ /- + oo 

D '(*' y) = mT) J 0 di ^ 8 ~ lT ( t ’^y) ■ ( 2 - 31 ) 

Similar results hold for Tr A _s ; for example, using the cylinder trace T(t) of Eq. 
(I2.29|) . we obtain 


p+oo 


Tr A~ a = 


r(2 s) 


dt t 28 - 1 T{ t) . 


(2.32) 


2.8 Analytic continuation of D s via complex integration. Let us assume 


T( t; x, y) = — J(t; x, y) , (2.33) 

for some q G Z, where the map J : [0, +oo) x x fl —> R admits an extension 
J : U{[ 0, Too)) xQxQ -)• C (W([0, Too)) C C is an open neigbourhood of [0, Too)); 
besides, for fixed x,y Gfi, assume the function t G 7/([0,Too)) (->• J(t;x,y) to be 
analytic and exponentially vanishing for 9?t —y Too. Then, one can infer from Eq. 
(12.31 j) that 

D.(x, y) = 1 TTLzgf) Jjte-'n t; x, y) (2.34) 

where denotes the Hankel contour, that is a simple path in the complex plane 
that starts in the upper half-plane near Too, encircles the origin counterclockwise 
and returns to Too in the lower half-plane (see Part 1 for an illustration of and 
for a precise definition of the complex powers of t). 

Eq. (12.341) gives the analytic continuation of D s to a meromorphic function of s 
on the whole complex plane, possibly with simple poles for s G {qj 2, (q — 1)/2, ( q — 
2)/2, ...} \ {0, —1/2, —1, —3/2,...}. For half-integer values of s, the integral in Eq. 
(I2.34p can be computed via the residue theorem; in this way for example, for s = 
—nj 2 and n G {0,1, 2,...}, we obtain 

D_ f (x,y) = (—l) n r(nTl) Res(t _ ^ n+1 ^ T(t; x, y); o) . (2.35) 


Using the modified cylinder kernel T , we deduce a variant of Eq. (12.341) : 

B.(x. y) = - e ' 2 ”" 2 r / 2 ~ 2s> l[ di i 2 " 2 n t; X, y ) . (2.36) 
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Again, for half-integer s the above analytic continuation can be computed explicitly 
by the residue theorem; more precisely, for n G {—1, 0,1, 2,Eq. (I2.36P gives 


D- f (x,y) = (-ir +1 r(n+2)Res(t-(" +2 )T(t;x,y);0) . (2.37) 

Similar results hold for the spatial derivatives of D s ; moreover an analogous dis¬ 
cussion can be made for the trace Tr*4“ s starting from Eq. (j2.32j) and using the 
cylinder trace T(t). Assuming the latter to admit a meromorphic extension to a 
neighborhood of [0, +oo), with only a pole at t = 0 and vanishing exponentially for 
9ftt —» +oo, we have 


TvA~ s = C 2s ) f di t 2s - l T(t) . (2.38) 

2m J a 

For s = —n/2 and nG {0,1, 2,...}, the above relation and the residue theorem give 
Tr A n/2 = (—l)”r(n + l)Res(V (ri+1) r(t);0) . (2.39) 

2.9 The case of product domains. Factorization of the heat kernel. Let 

A := — A + V and consider the case where 

0 = fi!x0 2 3x = (xi, x 2 ), y = (yi, y 2 ) , (2.40) 

y(x) = y 1 (x 1 ) + y 2 (x 2 ) ( 2 . 41 ) 

(fi a c R da is an open subset, for a G {1,2}; d x + d 2 = d); assume the boundary 
conditions on dfl to arise from suitable boundary conditions prescribed separately 
on <9fR and dfl 2 so that, for a = 1,2, the operators 

Aa := —A 0 + v(x a ) (2.42) 

(with A a the Laplacian on f2 a ) are selfadjoint and strictly positive (or at least, non¬ 
negative) in L 2 (Q a ). Then, the Hilbert space L 2 (Q) and the operator A can be 
represented as 

L 2 (fl) = L 2 (fA) <g> L 2 (fl 2 ) , A = Ai®1 + 1®A 2 . (2.43) 

This implies, amongst else, that the heat kernels K( t; x, y) := e t4 (x, y), K a ( t; x 0 , y a ) 
e tAa (x a ,y a ) (a = 1,2) are related by 

A'(t;x,y) = A'i(t;x 1 ,yi)/L 2 (t;xi,x 2 ) . (2.44) 
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2.10 The subcase of a slab: reduction to a lower-dimensional problem. 

By definition, we have a slab if 

x K d2 3 x = (xi, x 2 ), y — (yi, y 2 ) , V(x) = V (x x ) (2.45) 


with O i C H dl an open subset (di+d 2 = d), and if the boundary conditions prescribed 
for the field only refer to dfli x R d2 . We write D s (x, y) for the Dirichlet kernel of 

A := —A + V(xi) at x = (xi,x 2 ), y = (yi,y 2 ) G ; -Ds^(xi,yi) is the Dirichlet 
kernel of A\ := —A 1 + D(x 1 ). There hold 


Du± i(x,y) 


r( 


u~d2zkl ' 

2 


y=x 


(4vr)-i2/ 2 r( 


mb 1 ' 

2 ' 


D 


(i) 


u — ^2 A1 
2 


(xi,yi) 


yi=x! 


(2.46) 


d ^i D H i(x ’ y ) 


= o . P 


j LJ u-\-1 


y) 


y=x 


y=x 


= d i jDu+i (x. y) 

VaVb — V 


= 0 


y=x 


for (a, b ) = (1, 2) or (a, b) = (2,1) and i G {1, ...d a }, j G {1, ..., d b }; 

d-2 

d ziw iDu± i(x,y) - - 


T(^f±i) 


y=x 


(47r) d 2/2 r(2±i) ^ 


9^i>£.U(xi,yi) 


yi=xi 


U —<^2 — 1 


for z,w E {x, y} and i, j G {1,..., di} ; 

2 ) 


<1; v jD s ±i(x,y) = - d < jDu+i(x,y) = - d 4 ,D«+i (x, y) 

x 2^2 2 y=x X 2 X 2 2 y=x 2 


^ (47r) d2 / 2 2T(^-) ^ > u ~ d ^ 1 ( Xl ’^ 


yi=xi 


y=x 

for i,j G {1,d 2 } . 


The above relations, along with Eq.s (I2.16ti2.18ll . imply 

(0|7)“(x)|0) = 0 for i,j G {di + 1, ..., d}, i^j; 
(0|f)“(x)|0) = (0|fj)(x)|0)=0 foriG{l,...,d 1 },iG{d 1 + l,...,d} . 

Clearly, analogous relations hold for the renormalized VEV (0|T JJ (x)|0) r . eri 


(2.47) 


(2.48) 


(2.49) 


(2.50) 


2.11 Reduced energy for a slab configuration. Consider the slab configura¬ 
tion of subsection 12.101 the reduced regularized energy (i.e., the total energy per 
unit volume in the “free” dimensions) is 


£“:=/ dx 1 (0|fS,|0) = ^ + 5?. (2-51) 

Jn i 

The second equality is proved after defining the regularized reduced bulk and bound¬ 
ary energies, which are 


E\ L := 


« u r( 

2(4vr 


p/ u—d 2 — 1 \ /> 

F#y)l dXlD ^ (Xl ’ Xl) 


k u r( 


2 

XL — d>2 — 1 

2 


2(47r) d2 / 2 T(^i) 


3TT Tr -Ai 


<^2 + 1 — u 


(2.52) 
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(2.53) 


. U -p/ Ti-d2 + l \ /I \ C r) 

B ' = (Spkxisj (i ~vLf aM m D ^ uyi} 

One has = 0 for Ox bounded and either Dirichlet or Neummann boundary 
conditions on dfli. Assuming the functions (I2.52p (j2.53f) to be finite and analytic for 
suitable w G C, we define the renormalized, reduced bulk energy by the generalized 
(or restricted) zeta approach: 


E\ en := RP 


E 

u =0 


u 

1 


(or K" : " - H\‘ 



(2.54) 


2.12 Pressure on the boundary. This is the force per unit area produced by 
the quantized field inside h at a point x G <9fl We first consider, for large, the 
regularized pressure p“(x) with components 


Pi( x ) ■= (0|^(x)|0)rP(x) ; (2.55) 

here and in the remainder of this paper, n(x) = (n*(x)) is the unit outer normal at 
x G dfl. For Dirichlet boundary conditions the above definition implies 


p u M) = * u 


\ S U ° X ‘ 0 y‘ + l d x V 


D u+i(x,y) 


n J (x) 


y=x 


We can define the renormalized pressure by analytic continuation as 


(2.56) 


P r(x) := RP 


u =0 


pm 


(2.57) 


(first compute the regularized pressure at the boundary, and then analytically con¬ 
tinue at u = 0); alternatively, we could put 


iT(x) 


lim m 

x'eS!,x'^x 




Cx / )|0) r 


n J x 


(2.58) 


(first renormalize the stress-energy VEV at inner points of D, and then move to the 
boundary). 

Prescriptions (12.57P (12.58P do not always agree (for a counterexample, see Section 
[5]). In Part I we conjectured that the two approaches agree when both of them give 
a finite result (this is true in all the examples of the present Part II). 

2.13 The Hilbert space when 0 is an isolated point of cr(^4.); the case of 
Neumann and periodic boundary conditions. Assume A = —A+P, acting on 
L 2 (f2), to have its spectrum contained in [0, +oo), with 0 an isolated point (hence a 
proper eigenvalue); in this case, we replace the basic Hilbert space L 2 {Vt) with 

L 2 (D) := (ker A) 1 ' ( C L 2 (D)) . (2.59) 
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The restriction of A to Lg(12) is selfadjoint and strictly positive; we take Lg(12) as 
the basic space even for the field quantization. 

For example, if A = — A, 12 is bounded and either Neumann or periodic boundary 
conditions are prescribed on <912, one finds 


Lim 




dx/(x) 



(2.60) 


For slab configurations where 12 = 12x x R d2 and Neumann or periodic boundary 
conditions are prescribed on <912i x R^ 2 , we set {A\ is the reduced operator in L 2 (12i)) 


L&Sh) ■■= (ker Ai) 2 - 


/ e T 2 (n 1 ) 



dx! /(Xi) 



the basic Hilber space for the full theory on 12 is Lq(12i) <g) L 2 (H d2 ). 


(2.61) 


2.14 The case where 0 is in the continuous spectrum of A. Assume the 
fundamental operator A = —A + V to be non-negative (cr(A) C [0, +oo)), with 0 in 
the continuous spectrum of A. The approach we consider in this case is to represent 
A as 

A := “ lim ”A £ (2.62) 

where, for £ G (0, £o), A e is a selfadjoint operator in L 2 (12) with cr(A s ) C [e 2 , +oo) . 
We define a deformed, smeared field operator (j) £U := ( y n~ 2 A e )~ u ^(l) and a deformed, 
regularized stress-energy tensor operator T £ ™, whose VEV is 

(Ofe0r)|0> = 


— ^ (dxiiyv + d x v yl f) — ^ — — 2(j r/ /t „ (d x d y \ + V^j — £ (d x n x v + d y ^ y f) 
similarly, for x G <912, we consider the deformed, regularized pressure 


PT(x) : =<0|^(x)|0)n 2 (x 


In the end, we put 


{0\T^(x)\0) ren : = lim RP (0|T™(x)|0) ; 


£—>0 + 


/F en (x) := lim RP 
1 £^0+ 


u=0 


u =0 


p 


f(x) 


y=x 


(2.63) 

(2.64) 

(2.65) 

( 2 . 66 ) 
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We could alternatively put 


K en ( x ) : = ( , Jjm <0|7; i (x , )|0) ren > ) ra J (x) (2.67) 

\ x'Ga2,x'—>X / 

((OlTj^x^lO)^ is defined via Eq. (I2.65P ). For the VEV (I2.63P we have an expression 
analogous to (12.16112.181) in terms of the deformed Dirichlet kernel 


£h(x,y) := A s e (x, y) = (<5 x |^<J y ) , 


( 2 . 68 ) 


with s = (u ± l)/2. For the renormalized stress-energy VEV (12.65P we have an 
expression of the form (I2.211l2.23p . where 



lim RP 

£—> 0 + 



lim RP 

e—►()+ 


s=± 2 


(fi; 2sTl LL(x,y)) , 


(2.69) 


and analogous definitions hold for the spatial derivatives in the cited equations. 

In Part I we showed that two useful choices of A e are the following (K £ , T £ and K, T 
denote the heat and cylinder kernels associated to A e and A, respectively): 


A e :=A + e 2 => Jl e (t;x,y) = e £2t Jl(t;x,y) ; (2.70) 

A £ \= (\fA + e) 2 => T £ {i ; x, y) = e~ £i T(t; x, y) . (2.71) 

Let n G {—1, 0,1, 2,...}. If the modified cylinder kernel T(t;x, y) of A admits a 
meromorphic extension in t to a neighborhood of [0, +oo) fulfilling 

|T(t;x, y)| ^ C |t| _a_Tl+1 for 5Rt —> +oo and some C, a >0 , (2.72) 


then, assuming A e to be as in Eq. (I2.7ip . we have 


y) := lim RP 

2 £-> 0 + 


a (V+"CJ(x, y)) = 


(_ 1 )”+ 1 r(n+2) Res(r (n+2) f(t; x, y); o) . 


(2.73) 


With analogous hypoteses for d zw f(t; x, y) (z,w any pair of spatial variables), sim¬ 
ilar relations for the spatial derivatives can be deduced. 

2.15 Some variations involving the spatial domain. Configurations involving 
periodic boundary conditions are formulated rigorously describing fl in terms of tori; 
e.g., the domain 12 = (0, a) d with periodic boundary conditions is viewed as the torus 

Tf := K d /(aZ) d ~ (R/aZ) d 0). 

3 The considerations of subsection 12.131 for the periodic case are easily rephrased in terms of the 
torus T d a (see footnote 21 in Part I). 
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Sometimes (see Section EH) we employ on some set of curvilinear coordinates 
(q l )i=i,..., d = q, inducing a set of spacetime coordinates q = (g 4 *) = (t, q); the 
spatial and space-time line elements are, respectively 


d£ 2 = aij(q)dq l dq :j ; ds 2 = —dt 2 + d£ 2 = g^ v {q) dq^dq u , 

9oo ■= -1 , gio = goi ■= o , Sy(g) := aij( q) for i,j G {1, ...,d} . 

The analogue of Eq. (12. 8 p in the coordinate system (g M ) is 

f; u := + -2£0“o (2.75) 

(V M is the covariant derivative induced by the metric (I2.74p ). For any scalar function 
/ there hold (yh are the Christoffel symbols for the spatial metric (ay(q))) 

V,/ = d,f , V y / = Dijf = duf - 7£A/ , f2 76) 

v 0i / = 9o(^/) = ^(a 0 /) = Vj 0 / , Voo/ = 9oo/ • 

Many results in the previous subsections are readily adapted to the variations con¬ 
sidered in this subsection for the space domain. 


3 The case of a massless field between parallel 
hyperplanes 

3.1 Introducing the problem for arbitrary boundary conditions. As men¬ 
tioned in Section 6 of Part I, the segment configuration can be considered as the 
d — 1 case of a general, d-dimensional configuration with two parallel hyperplanes; 
this is the subject we are now going to analyze (with no external potential). So, we 
assume 

n := (0, a) x R d l , a > 0 , V = 0 ; (3.1) 

these choices correspond to a massless scalar field confined between the two parallel 
hyperplanes ( 0 ) 

7r 0 = {x e R f/ | x 1 = 0} , TT a = {x G R f/ | x 1 = a} . (3.2) 

4 There holds comments analogous to the ones in footnote 18 on page 49 of Part I; namely, 
in place of the standard Cartesian coordinates x = (ad), we could have used the set of rescaled 
coordinates (best fitting the features of the present configuration) 

x* = (xl )i=i,...,d , with x\:=x 1 /a, x\ := ad for i G {2,..., d} . 

Also in this case, we choose not to employ the above rescaled coordinate system in order to render 
comparison with kwnown results easier. 
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As we did in Section 6 of Part I for the segment configuration, we are going to 
consider separately the cases where the field fulfills Dirichlet, Neumann or periodic 
boundary conditions on the hyperplanes 7r 0 ,7r a . Throughout this section we assume 


(3.3) 


d odd , d ^ 3 ; 


this hypothesis is purely technical and will be motivated later (see the comments 
after Eq. (13.61) ): note that the case d — 1, here excluded, has been already discussed 
in Part I. 

In passing let us notice that, for d — 3 and Dirichlet boundary conditions (see 
subsection 13.6[) , the above configuration is the one most typically considered when 
dealing with the (scalar) Casimir effect 0125] 130]. The case with Dirichlet boundary 
conditions on one plane and Neumann conditions on the other (discussed in sub¬ 
section 13.71) was originally considered in the electromagnetic case by Boyer |5j, who 
derived the total energy; later, computations of the total energy and boundary forces 
for a (massless or massive) scalar field at both zero and non-zero temperature were 
performed by Pinto et al. [ID) [28] and Santos et al. [3_2j (see also Ban. Finally, 
let us also mention the monography by Fulling [22] where the stress-energy VEV for 
the model with periodic boundary conditions is given; see also nans] and, again, 
0121 for the derivation of the total energy in the same configuration. 

For any one of the boundary conditions mentioned above, keeping in mind the 
considerations of subsection 12.131 we proceed in the manner explained in the fol¬ 
lowing subsection. Before moving on, let us remark that, due to the results on 
slab configurations reported in subsection 12.101 we just have to study the reduced 
one-dimensional problem based on 


- - i • (0, u) C R. , wT-i . d x i x i , 


(3.4) 


keeping into account the boundary conditions in x l = 0 and x 1 = a ; in consequence 
of this, we can resort to the results of Section 6 in Part I. 

In subsections 13.2113.51 we will present some general results on the configuration 
under analysis, holding for all the types of boundary conditions mentioned before. 
In subsections 13.6113.91 we will consider specific boundary conditions, with a special 
attention for the case d = 3. 

3.2 The reduced Dirichlet and cylinder kernels. According to Eq.s (12.461 
I2.49p . the basic ingredients for the analysis of the d- dimensional problem are the 
Dirichlet kernel Ds of the reduced 1-dimensional problem at s = (u — d)/ 2, and 
its derivatives at s = {u — d + 2)/2. On the other hand, these functions can be 
expressed in terms of the 1-dimensional cylinder kernel Td)(t; x 1 , y 1 ), which has been 
determined in Section 6 of Part I (where it was indicated simply with T(t; x 1 , y 1 ))', 
let us recall that this kernel and all its derivatives, when evaluated on the diagonal 
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y 1 = x 1 , are meromorphic functions of t in a neighborhood of the positive real half¬ 
axis with a unique singularity in t = 0, and they vanish exponentially for 9?t —y +oo. 
For the reduced Dirichlet kernel and for the derivatives in which we are interested, 
Eq. (12.34)1 yields the following expressions: 


D 


(i) 

u — d 
2 


/ 1 1\ e d) T{d+l U ) f u-d-l rp(l)(, 1 1\ (O r\ 

(x ,y ) = - — - dit T (, (i;x,y)] (3.5) 


dzwPu-d+2 y 1 ) = 


l u e- i < u - d ^T{d-l-u) 


2vr i J 

for z,w € {x 1 , y 1 } 




(3.6) 


Consider the above relations along with Eq.s (12.46112.49]) . relating the 1-dimensional 
functions to the d- dimensional Dirichlet kernel Du- 1 (x, v) or to the derivatives of 
D «+i(x, y); these equations yield the sought-for meromorphic continuations in u 
to the whole complex plane. By direct inspection of the expressions thus obtained 
it appears that, with the assumption (13.31) on d, these meromorphic continuations 
are analytic for u in a neighborhood of the origin; so, the zeta strategy for renor¬ 
malization is implemented by simply setting u — 0. In conclusion, we have the 
following integral representations for the renormalized Dirichlet kernel and for its 
renormalized derivatives: 


£>i(x,y) 


c d r dt 


y=x 


27 ri ./fj t d+1 


T (1) (t;x 1 ,7/ 1 ) 


(3.7) 


d zw D i(x, y) 


C, 


dt 


y=x 


(d— l)2m t d 


r-rd zw T { - l) ( tjx 1 ,?/ 1 ) for z,w G {x 1 , y 1 } ;(3.8) 


diiD i(x,y) =-dijD i(x,y) = - d yi jDi (x, y) 

X 2V2 2 y_ x 2 y— x V2V2 2 


y— x 


= 5,: 


Cb /• dt 


* 27T7 ^ D+ 

for the sake of brevity, in the above we have set ^ 


Y T^^(t; x 1 , t/ 1 ) for i, j G {1,..., d —1} ; 


(3.9) 


<?«*:= 


(-7T)' 


d~l 

— r 



(3.10) 


Eq.s (13.7H3.9p are completed with Eq. (I2.47p . stating the vanishing of certain mixed 
derivatives. 

Finally, recall that Eq. (j2.52[) (here employed with d 2 = d—1) for the regularized, 
reduced bulk energy (see the subsequent subsection 13.41) requires the evaluation of 

5 To justify the expression (13. 101) . some identities regarding the gamma function must be used. 
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the trace TrA 1 2 . To this purpose, we first consider the one-dimensional cylinder 
trace T^)(t), which has also been determined in Section 6 of Part I (where it was 
indicated simply with T(t)); we recall that, for all the boundary conditions con¬ 
sidered in the following applications, the map t T^(t) admits a meromorphic 
extension to a neighborhood of [0, +oo) which only has a pole at t = 0 and vanishes 
exponentially for JR —* +oo. 

A discussion similar to the one carried over above for the Dirichlet kernel allows us 

d—u 

to derive an explicit expression for the analytic continuation of Tr^l 1 2 at u = 0; 
more precisely, using for the reduced problem a relation analogous to one in Eq. 
(I2.38p . we conclude 

TvA d/2 = {-l) d T(d+l) r dii -( d+ i) T (i)( t ) _ (3.H) 

2m 


3.3 The stress-energy tensor. Substituting the relations (12.47P and (l3.7H3.9p 
into Eq.s (l2.16H2.lgp . we straightforwardly deduce the contour integral represen¬ 
tations for the non-vanishing components of the renormalized stress-energy VEV; 
moreover, due to the meromorphic nature of the cylinder kernel (and of its deriva¬ 
tives), the resulting integrals along the Hankel contour can be explicitly evaluated 
via the residue theorem. The final expressions for the (non-zero) components of the 
renormalized stress-energy VEV are the following ones: 


(0|f 00 (x)|0) ren = -C d Res ( r (d+1) 


d -2 


4 d 


€ - TT )dT ( ' 1 \i-,x\y 1 ) + 


(l <9 xV T (1) (t;xV) 


d -1 V 4 


(0|Tn(x)|0) ren = —C d Res Q - ^jd T (1) (t; x 1 , y 1 ) + 

+ Q d x i y i - ^ X i 3 .i'jr (i) (t;x 1 ,y 1 ) 


;0 ; 


(0|^(x)|0) re7 , = 8 tj C d Res ( r (d+1) 


e- 


d-2 


4 d 


+ JZI (\ -t) d xiy'T {1) { t^x^y 1 ) 


(3.12) 


;0 ; 


(3.13) 


d T" (1) (t; a? 1 , y 1 ) + 

; 0 J for i,j e {2,..., d} . 


(3.14) 


We repeat that, in the above, d is an arbitrary odd dimension > 1. Starting from 
subsection 13.61 for each one of the previously mentioned boundary conditions we will 
report the explicit expressions for the the stress-energy components arising from Eq.s 
(13.12113.14P in the case of spatial dimension d = 3; again, we will give the final results 
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in the form described in subsection 12.31 (see, in particular, Eq. (I2.13p ). noting that 
(12.14P gives 

£3 = \ • (3.15) 

6 

Some additional details related to these computations will be given, as examples, in 
the case of Dirichlet and periodic boundary conditions. 

3.4 The reduced energy. Let us recall again that this is the energy per unit 
volume in the free dimensions and that it can be expressed in terms of the reduced 
bulk and boundary energies (see subsection 12. lip . 

The general identity (I2.52p allows us to represent the reduced bulk energy in terms 
of the renormalized trace Tr , corresponding to the reduced operator on the 
other hand, Eq. (13. lip gives an explicit expression for the latter quantity in terms of 
the reduced cylinder trace Td'(t). Evaluating the integral along the Hankel contour 
in the cited equation via the residue theorem (compare with Eq. (12.391) 1. we conclude 

E T n = Re s (r( d+ 1 >r<l>(t); 0 ) . (3.16) 


Let us also mention that, for any one of the boundary conditions to be considered in 
the following, the (regularized and renormalized) reduced boundary energy always 
vanishes identically. 

Finally we point out that, at least in spatial dimension d — 3, the results obtained via 
Eq. (I3.16P for the renormalized, reduced bulk energy coincide with the integral over 
the interval (0, a) of the conformal part of the corresponding renormalized energy 
density (0|T 0 o|0) ren ; on the contrary, the non-conformal part of the latter appears 
to diverge in a non-integrable manner near the planes 7r 0 ,7r a . These facts closely 
resemble the ones pointed out for the total energy of the segment configuration in 
Section 6 of Part I; as in that case, they will be checked by direct computation in 
subsections I3.6H3.9I 

3.5 The boundary forces. The situation we meet in the present situation is of 
the kind described in general in subsection 12.121 and already faced for the segment 
configuration in Part I: in principle, we can define the renormalized pressure on the 
plates 7r 0 ,7 T a in two alternative ways. 

Let n(x) denote the unit “outer normal” at points on the boundary, so that n(x) = 
(—1, 0,..., 0) on 7r 0 and n(x) = (1, 0,..., 0) on 7r a . Then, on the one hand we can put 


K en (x) := <0|7™(x)|0> 


n-Tx' 


u =0 


= Sn (0|T“i(x)|0) 


u =0 


(3.17) 
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on the other hand, we have the alternative definition 


PT(x) := , Jim (0|7ij(x / )|0) ren n J (x) 

' x'Gi2,x '—• 


= M , 1™ (°l T ll( X, )|0)r 

.x'Gr^x'—>x 


(3.18) 


As a matter of fact, definitions <13.171) and <13.181) will be found by explicit com¬ 
putations to yield the same result for any one of the boundary conditions to be 
considered in the next subsections. 

3.6 Dirichlet boundary conditions. Let us first consider the case where the 
field fulfills Dirichlet boundary conditions on both the hyperplanes 7r 0 ,7r a , meaning 
that 

0(i, x) = 0 for t G R, x G 7r 0 or x G 7r a . (3.19) 

Recall that in this case the cylinder kernel associated to the reduced problem is (see 
Eq. (6.20) in Part I) 


T (1) (t;x 1 ,r/ 1 ) 


1 cos(^(x 1 — y 1 )) — e a 1 

2a cosh(^t) — cos(^(x 1 — y 1 )) 


cos(-(x 1 +r/ 1 )) — e 
cosh(yt) —cos(^(:r 1 +r/ 1 )) 


(3.20) 


To obtain the renormalized stress-energy VEV in any spatial dimension d, it suffices 
to substitute Eq. (13.20)1 into Eq.s (13.12113.14)) and to compute the residues therein. 
Let us explicitate the final results for d — 3; in this case the residues in (13.12113.14)) . 
involving T^ii^x 1 ,y l )\ y i =x i, can be derived from the t —* 0 expansion 


T (1) (t; x 1 , y 1 ) = 

y 1 =X 1 

1 _ irp-siipfx 1 )) 7 t 3 (15(2 + cqs(^ l x 1 )) — sin 4 )^ 1 )) < 3 ' 21 ) 

7rt + 12a 2 sin 2 (- x 1 ) 720a 4 sin 4 (£ x 1 ) 1 

Proceeding in a similar manner where the spatial derivatives of T' 1 ) appear, we 
obtain this final result for the d = 3 renormalized VEV of the stress-energy tensor: 


(0|f^(x)|0) ren 


H,is= 0,1,2,3 


/-I 0 0 0\ 

0-300 
0 0 10 
\ 0 0 01/ 



/-1 0 0 0\ 
0 0 0 0 
0 0 10 
\ 0 001/ 




A 


1440a 4 ’ 


B^x 1 ) 


7r 2 3 — 2 sin 2 (^ x 1 ) 
sin 4 (- x 1 ) 

'a / 


for x 1 G (0, a) . 


(3.22) 


This is a classical result, whose earlier derivations used point splitting methods 
rather than zeta regularization (see, e.g., the paper by Esposito et al. [16] or the 


21 


























already cited monographies [211221 [25]) ; in our previous work ra the expression 
(I3.22p was obtained via zeta regularization, but the required analytic continuations 
were derived by ad hoc considerations, strictly related to the peculiar conhguration 
under analysis. On the contrary, here the analytic continuation arises automatically 
from the general schemes that we have developed for an arbitrary geometry (0). 

Next, let us pass the reduced bulk energy. To this purpose, recall that the cylinder 
trace associated to the reduced problem under analysis is (see Eq. (6.25) in Part I) 

T (I) (t) = —; (3.23) 

e« — 1 


now, the general relation (13. 1611 allows us to infer, for d = 3, 


ren 


1440a 3 ' 


(3.24) 


Finally, using either definition (13. 17j) or (13.18j) . we obtain for the non-vanishing 
component of the pressure on the planes hq and 7r a the following expressions, re¬ 
spectively: 


pr(x) = 3kL , pj"*(x) 

7TQ 


= -3A 


(3.25) 


where A as in Eq. (13.221) . This means that in the present case, the forces on the 
boundary planes produced by the held in the interior region are attractive (0). 


3.7 Dirichlet-Neumann boundary conditions. Consider now the parallel hy¬ 
perplanes conhguration where the held fulfills Dirichlet and Neumann boundary 
conditions, respectively, on the hyperplanes 7r 0 and 7r a ; explicitly, 


</>(£, x) = 0 for (t,x) e R x 7r 0 , d x i 0(i, x) = 0 for (t,x) e Rx7r a . (3.26) 


The cylinder kernel associated to the reduced operator A\ is (see Eq. (6.30) in Part 

I) 


T {1) ( tjarSj/ 1 ) 


1 r sinh (^t) 

a cosh(^ t)— cos(^(x 1 — y 1 )) 


sinh(^t) cos(^(o; 1 +|/ 1 )) ~ 
cosh(f t)-cos(^(x 1 +7/ 1 ))_ 


(3.27) 


6 In connection with the present approach, see also [15] . 

7 Clearly, we are not taking into account any effect related to the outer region (see the comments 
at the end of subsection 12.121 of Part I). As a matter of fact, the forces produced by a massless 
scalar field in this region, fulfilling either Dirichlet or Neumann boundary conditions on the planes, 
vanish identically for d = 3; this result can be derived using the methods that will be developed in 
the following Section |4] 
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Employing this kernel along with relations (13.12113.lTp . we can evaluate the renor¬ 
malized VEV of the stress-energy tensor; in particular, for d = 3 we obtain 


(0|T^(x)|0) r 


= A 


0 , 1 , 2,3 


A = 


7n 2 


11520a 4 ’ 


^(a; 1 ) = 


/I 0 0 0 \ 

0 3 0 0 

0 0-10 

\o o 0 - 1 / 

7T 2 23cos(^£ 1 )+cos(^ L ir 1 ) 


-u-gjsc* 1 ) 


(-1 0 0 0 \ 
0 0 0 0 
0 0 10 
\ 0 001/ 


64 a 4 


sin 4 (- x 1 1 


for x 1 E (0, a) . (3.28) 


To the best of our knowledge, the full stress-energy VEV for the present configuration 
has never been given in the preceeding literature; only the global results on energy 
and pressure discussed hereafter have previously been considered (see the citations 
below). 

Let us recall that the reduced cylinder trace in this case is (see Eq. (6.33) in Part I) 


T (1) ( t) 


— t 

e2a 



(3.29) 


The above expression, along with prescription (13.16|) . allows us to determine the 
renormalized, reduced bulk energy; for example, for d = 3 we obtain 


jjjren 


7n 2 

11520a 3 ' 


(3.30) 


Concerning the non-vanishing component of the boundary pressure, both definitions 
(I3.17P (I3.18P give (with A as in Eq. (I3.28P ) 


K en (x) 


= -3A , 


7TQ 



(3.31) 


let us stress that, similarly to the results found for the segment configuration in Part 
1 (see Eq.s (6.27) (6.35) therein), the expressions in Eq. (I3.3ip have the opposite 
sign with respect to those in Eq. (I3.25p . holding in the case of Dirichlet conditions 
on both the planes 7r 0 and 7i a . This means that in the present Dirichlet-Neumann 
case, the forces on the boundary planes are repulsive (jf|). 

The results in Eq.s (I3.30P (I3.3ip agree with those reported, e.g., in [10], [28j and 
Santos et ah [32J (see also Ha). 

3.8 Neumann boundary conditions. Assume that 


d x i4>(t, x) = 0 for t e R, and x e 7r 0 or x e n a ■ 
8 Recall the considerations of the previous footnote [3 


(3.32) 
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Recall that in this case, according to the considerations of subsection 12.131 the 
reduced operator A\ must be viewed as acting in the Hilbert space Ljj( 0, a) of square 
integrable functions on (0, a) with mean zero (see Eq. (j2.60jl ). The cylinder kernel 
of A\ is (see Eq. (6.38) in Part I) 


T (1) (t; x 1 , y 1 ) 


J_ f cos(^(x—y)) -e * cos(f (x+y)) - e 1 

2 a [cosht — cos(^(x—y)) cosht — cos(^(x+y)) 


(3.33) 


Using this kernel along with relations (13. 12113.14)) . we obtain the renormalized stress- 
energy VEV; in particular, for d = 3 the residue evaluation yields 


(0|T MI y(x)|0) ren 


/i,^=0,l,2,3 



(-1 0 

0 

°\ 


(-1 

0 

0 

°\ 


0 -3 

0 

0 

A 1\ 

0 

0 

0 

0 

A 

0 0 

1 

0 

+ ({- 5 )sa) 

0 

0 

1 

0 


O 

O 

0 

V 


V o 

0 

0 

V 


(3.34) 


where A and B^x 1 ) are defined as in Eq. (13.221) . 

Let us pass to the reduced bulk energy; in subsection 6.8 of Part I we noticed 
that in this case the cylinder trace T^(t) coincides with the one of Eq. (13.231) . 
corresponding to the case of Dirichlet boundary conditions. In consequence of this 
the renormalized, reduced bulk energy for d — 3 is the same as in the case of Dirichlet 
boundary conditions (see Eq. (13.24ft 1: 


1 1440a 3 ' 

Furthermore, concerning the pressure on the boundary, both definitions (I3.17|) and 
(I3.18p also give the same result we derived in the case of Dirichlet boundary condi¬ 
tions (see Eq. (I3.25[) : again, A is the coefficient of Eq. (13.22[) 1 


pF(? 


= 3H 


K en (x) 


7TQ 


= -3A . 


3.9 Periodic boundary conditions. The last case we consider is the one where 

4>{t, 0,x 2 ,... ,x d ) = 4>(t,a,x 2 , ...,x d ) , d x i4>(t, 0, x 2 ,... ,x d ) = d x i4>(t,a,x 2 , ... ,x d ) 

for t, x 2 ,..., x d e R . (3.35) 

Similarly to what was said for the segment with periodic boundary conditions, note 
that, as explained in subsection 12.151 this configuration would be more properly 
formulated in terms of a free scalar field on the flat manifold : = T^x R d_1 , 
where the first factor is the torus := R/(aZ). 
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As in the case of Neumann boundary conditions, the basic Hilbert space for the 
reduced problem is L^T 4 ) (see subsection 12.131 Eq. (I2.60p ). We know that the 
cylinder kernel associated to Ai in this case is (see Eq. (6.44) in Part I) 


T (1) (t; a; 1 , y 1 ) = 


cos (^(x 1 -? 1 ))-^* 
a [cosh(^t) - cos^x 1 -?/ 1 ))] 


(3.36) 


Again, we can employ Eq.s (13. 12113.141) to evaluate the renormalized stress-energy 
VEV. Differently from the previous subcases, this time the expressions appearing in 
intermediate steps of the required calculations are simple enough to be reported; as 
an example, let us focus on the evaluation of the component (0|T 0 o(x)|0) ren . First 
of all, note that 


T (1 /t; x 1 , x 1 ) = - 


a l 


coth ( — t 


d x \yiT^\i] x 1 , x 1 ) = —— coth (—t) csch 2 f—t 


So, after some simple algebraic manipulations, Eq. (I3.12p yields 


(3.37) 

(3.38) 


(0|?oo(x)|0) ren 


-—Res 

4a 


2 —(1—4 flrf 



+ 


2(1-4Q 
(d— 1) t d+1 



coth(^t)csch 2 (^t) 



(3.39) 


The function in the above expression, whose residue in t = 0 is required, is easily 
seen to be meromorphic with a pole of order d + 2 in t = 0; more precisely, its 
Laurent expansion is 


d(d-3)+4(d 2 -d-2)£ 1 (2-d)+4d£ 1 

47r(d—1) t d+2 + 4 a t d+1 

7r((2-d)+4dQ _1 _ 7r 3 (-(d 2 -3d-4)+4(d 2 -d-6)Q 1 

12 a 2 t d 180(d—1) a 4 i d ~ 2 1 ’ 


For example, for d = 3 Eq. (13.39j) yields (A 3 = — 1 / 7 r, see Eq. (13.10p ) 


(0|f 0 o(x)|0) ren = 


90 a 4 ’ 


(3.41) 


proceeding similarly for the other components of the renormalized stress-energy 
VEV (for d — 3), we obtain 


(0|T M „(x)|0) ren 


H,v= 0,1,2,3 


/-I 0 0 0\ 

7r 2 0 -3 0 0 

90a 4 0 0 1 0 

\ 0 0 01 / 


(3.42) 
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Notice that the renormalized stress-energy tensor 03.421) does not depend explicitly 
on the periodic coordinate x 1 . This comes as no surprise; indeed, it reflects the 
invariance of the theory under translations x 1 h->- x l + a (for arbitrary a G R). 

To conclude, let us consider the total bulk energy. First, recall that the reduced 
cylinder trace in this case is (see Eq. (6.47) in Part I) 

m = nr—r ; ( 3 - 43 ) 

e a — \ 

so, using once more the prescription 03.16p we obtain 

7T 2 

E ren = -- . (3.44) 

90a 3 V ; 

In conclusion, let us stress that the above expression for the total energy agrees with 
the results in 000 [25]. 


4 The case of massive field constrained by per¬ 
pendicular hyperplanes 

4.1 Introducing the problem. We consider a scalar field of nonzero mass m 
fulfilling either Dirichlet or Neumann boundary conditions on d\ orthogonal hyper¬ 
planes in d = d\ + g ?2 spatial dimension. More precisely, we are interested in the 
case 

12 = (0,Too) dl x R d2 , V(x) = m 2 (m > 0) . (4.1) 

The domain 12 is bounded by the hyperplanes {x 1 = 0},..., {x d = 0}; its boundary 
is the union of the faces 

7r n := {x G <912 | x n = 0} , (n G {1,..., di}) (4.2) 

and, for each one of them, either Dirichlet or Neumann boundary conditions are 
prescribed. As anticipated in the Introduction, the same configuration with at most 
three faces was considered by Actor et al. in OEj; when possible we will establish 
connections with these works, making direct comparison. 

Before proceeding, let us stress that also in this case we can use the results of 
subsection 12. 101 on slab configurations; because of this, we will consider the reduced 
problem based on 


12i = (0, +oo) dl , Ai := —Ai + m 2 , 
with the appropriate boundary conditions. 


(4.3) 
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4.2 The reduced heat kernel. In our approach, a basic step for the analysis 
of the reduced problem (14.31) is the computation of the heat kernel K^>(i ; xi, yi) 
associated to A\. The result is ||) 


K m ( t;xi,y,) 


e -m 2 i 

(47rt) dl / 2 





4t 


+ (y, n 6 


(*?■+W?) 

4t 



(4.4) 


where, for any n G {1,..., di}, a n G R is a parameter distinguishing between Dirich- 
let and Neumann boundary conditions on the face 7r n ; more precisely, 


J—1 for Dirichlet B.C. on 7i n 
{+1 for Neumann B.C. on 7r n 


Eq. il 1.11) can be re-written as 



(n G {1,..., d \}) . 


A' (1) (t;xi,yi) = 


e~ m2t 

(47rt)rfi/2 


£ 


77.—0 


1 

(di —n)\n\ 



(4.5) 


(4.6) 


9 Here is one way to derive Eq. Gd}. First, notice that a complete orthonormal system of 
(improper) eigenfunctions of Ai = — Ai + to 2 on f2i = (0, +oo) dl fulfilling the prescribed (either 
Dirichlet or Neumann) boundary conditions on cffh is given by 


1 di 

Fk(Xl) := ( 27r )^/2 II ( ^ e~ ik ^) , cu k := V|k| 2 +m 2 for k £ K, = (0, +oo) d i ; 

k ) n= 1 

here, for any n £ {l,...,di}, a n is defined according to Eq. (14.51) . Then, using the eigenfunction 
expansion (12.261) of the heat kernel, we obtain 


„ — m 2 t r+oo 

K m (i i XL yi) = n / <tfc<r “( 

„ —774 2 t dl «-f 


e ikx ? + e 


+Q!ne ifcyf j 


7 JJ / dfc e _tfe2 + an e <fc(*r+»D j ; 


( 2 7r) dl J 
v ’ n=l -y_ 


Eq. (14.41) follows by explicitly evaluating every single Gaussian integral in the product. 
10 This result depends on the identity 


(on+ b n ) — y) 


e n 


(d—n)!n! 

77=0 v ' cr€;Sd \l=l 



holding for any d G {1,2,3,...}, a n ,bn G R (n G {1, 2,..., d}), where by convention we intend 

rui a o-(») ;= rwn ^o-(j) ;= i • 
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Here and in the following S dl denotes the symmetric group with d\ elements and, by 
convention, the sums over i and j in (14. 6 p are zero for n = 0 and n = di, respectively; 
moreover, for any a G S dl , we put 

di 

: = &a( l ) for n G {0, , a duCT := 1 . (4.7) 

l=n +1 

4.3 The reduced Dirichlet kernel. In the following we will use Eq. (12.30(1 to 

express the Dirichlet kernel Z?i^(xi,yi) associated to Ai, along with its analytic 
continuation, in terms of K^\ Substituting the expression (14. 6 p for the heat kernel 
into Eq. (I2.30p . we obtain 


Di 1 ’ ( 


xi,yi) = 


di 

E 


(4 7 r)c* i /2 r(s) (di — n)\n\ 


E 

o-eSd 1 


^n,<T ^s,n,cr(^l 5 y 1J 


(4.8) 


where, for each nG{0, ...,di} and each aGS^, 


"+QO 


a w (xi, yi ) := / dti s -^- l e- mH -^^^ , 


(4.9) 


9l„ iCT (xi,yi) 


E 1 

. *=i 


x 


r W 


Vi 


(*)\2 


di 

+ E< 

j=n +1 


1/2 


y u) +i/E) 


(4.10) 


Clearly ^.^(x!, y x ) ^ 0 , so that convergence conditions for the integral in Eq. (14.9p 
can be readily infered; more precisely, if ^^(xi, yi) > 0 the integral converges for 
any s G C while if 9T ni(T (xi, yi) = 0 (which, in the interior of O, happens if and only 
if n — d i and yi = xi) it only converges for 


^ di 

s s > T 


(4.11) 


Under these assumptions for convergence, Eq. (14.9p is strictly related to a known 
integral representation of the modified Bessel function of the second kind K v (see, 
e.g., [27], pag.253, Eq.10.32.10); using this representation we obtain, for any x 1; y x G 
(0, +oo) dl , 


3i,n,a(xi, yi) = 2 2 1+1 s rn dl 2s 0 dl (m 2 C^ tr (xi,yi)) , (4.12) 

where, for the sake of brevity, we have put 

0, : [0, Too) -T C , 0„(z) := z v/2 K v {^) • (4-13) 
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Ill conclusion 


D * 1)(x ‘- yi) = (L)«X W E e yi)) • ( 4 - 14 ) 

v 7 w n=0 v ' o-eS dl 

The derivatives of Di 1 ' 1 of any order can be computed using the identity ((3) 


d n & ( 1' 

(z)={-- ) <5^ n (z) for nG {1,2,3,...}; 


(4.15) 


dz n v "' V 2 

moreover, in the cases where y^ = 0, we can resort to the relation 

0„(O) = 2 i " 1 r(z/) for V G C, 3*1/ > 0 . (4.16) 


4.4 The d-dimensional Dirichlet kernel. Using the previous results and Eq.s 
fl2.46H2.4Hl) . we obtain the following expressions for the Dirichlet kernel D s (x, y) of 
the d -dimensional problem (14.11) and for its derivatives, along the diagonal y = x: 


Du±i(x,y) 


2 2t1 2 u m dTl u 1 

y=x (27r) d / 2 r(“±l) 2 -^ (di-n)\n\ 


Y an,* g -ai (m 2 9 f l^ 0 .(xi, yj)) 

<*&S dl 


(4.17) 


yi=xi 


d z i w jDu± i(x,y) 


2 2 777, d 1 u dl 


y=x 




( 2,^/2 r(s±i) ^ (*-«)!«! 


o-eSd, 


(4.18) 


yi=xi 


for z,w E {x, y} and i,j G {1,..., d\} ; 


d x i v jDu± i(x,y) =-d i iDu+i(x,y) =-d i (x, y) 

x 2^2 2 y— x X 2 X 2 2 y— x * 2^2 2 


y~ x 




2 2 77 U + 1 u 


di 


*3 


unE rx—7(0 E“v ®*=f=i( m29l h( x i.yi)) 


(2^)4/2 r(!!±l) f- (d, -„)M 


(4.19) 


yi=xi 


2 7 71 =0 V 7 (leSdj 

for i,j G {1, ...,d 2 } • 

Note that, in each of the sums over n appearing in the above expressions, the terms 
corresponding to n G {0, l,...,di — 1} are analytic functions of u on the whole 


11 We already considered the map ©„ of Eq. (14.131) in Appendix D of Part I; therein we also 
showed how to derive the relations (14.151) (14.161) . 
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complex plane since 91); ^ (x x , x x ) > 0 for any x x G (0, +oo) dl . On the contrary, the 
terms corresponding to n = d\ deserve particular attention; indeed, ^(xi, x x ) = 0 
for any x x G (0, +oo) dl so that, in order to evaluate these contributions, we have to 
resort to Eq. (14.161) (also recalling Eq. (14.151) 1. In this way we obtain 



(4.20) 



(4.21) 


for z,w G {x, y} and i,j G {1, ..., d{\ . 

In principle, the above equations hold with the limitations on u arising from Eq. 
(14.161) : more precisely, Eq. (I4.20p holds for > d -Fl and Eq. (I4.2ip for > d+1. 
However, the right-hand sides of these equations are well defined and analytic on 
the whole complex plane with the exception of simple poles placed at 


u = d+ l-2£ for £ G {0,1,2,...} ; 


(4.22) 


this remark gives the meromorphic continuation in u of the functions in Eq.s (I4.20j) 
(14. 21 p and, consequently, of the terms with n = di in Eq.s (14.17114.191) . 

4.5 The stress-energy tensor. Using Eq.s (12.16112.181) along with Eq.s (14.171 
14.191 (and Eq.s (l4.20tl4.2TT) for the terms with n = di), we obtain the analytic 
continuation to a meromorphic function of each component of the regularized stress- 
energy VEV, required in order to implement the local zeta approach. In particular, 
we have the following expressions for the non-vanishing components: 



yi=xi 
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< 0 |ly(x)| 0 ) = ( 0 |T“(x)| 0 > = 


2 V m d+l nn\~ u 1 


Jl\ 


( 27r )^r(H±l) \ K J Z^( dl - n )\ n 


M±ii ( K. 1 E (d,-n.) In! ^ 


QV 




+ <g ^i (?n 2 ^n,a(xi,yi)) + g 2= d±i (m 2 9^ [ff (xi,yi))^ + 


m 


-2 


i -t) tijd^dyl + 9 xWi ~ e 0u-|+i (mX i<T (x 1 , yi )) 
for i,j e {1,..., di} ; 


yi=xi 

(4.24) 


(oK“W|o) = <o|T),?j(x)|o) = -s. 


f. 1 —u 

2 2 rri 






-'a? 


S 7 E 


(2x)^r(5±i)W £-(*-«)!«! 


E 


Q^r 


d — c?i — 1 — u 


n =0 
2(Y>2 


<*&S dl 


(d-di+l-u)£ g .-j-i (mX t r(x 1 'Yi)) + 


+ (t-£ (! + m d^dy^&u^+iim Tt^x^ yi )) 


yi=xi 


for i, j e {di + 1,..., cZi + c ?2 = d} . (4-25) 

The renormalized VEV of the stress-energy tensor is obtained sending u to zero in 
the above expressions; the only singularities appear in the terms corresponding to 
n — di, which must be treated resorting to Eq.s (14.20114.221) . 

The conclusions are the following: 

i) For d even each component of the regularized stress-energy VEV is an analytic 
function of u near u = 0; thus its renormalized version is obtained via the restricted 
zeta approach, i.e., by simply evaluating Eq.s 04.23H4.25]) at u — 0. 

ii) For d odd the regularized stress-energy VEV has a simple pole in u — 0, so that 
we have to resort to the extended zeta approach and consider the regular part at 
u — 0. 

The manipulations indicated in i) are trivial; the ones indicated in ii) could be 
performed in principle for an arbitrary odd dimension, but the final expressions are 
too lengthy to be reported here. For this reason, we prefer to exemplify ii) in two 
special cases with d — 3 (see subsections 14.81 and 14.91) . 

4.6 The boundary forces. As in the previous section, following the general 
framework of subsection 12.121 we can give two alternative definitions for the pressure 
acting on the boundary of the spatial domain . 

Let us consider a point x e <9f2; if d\ > 1 we exclude x to be on a corner, where the 
outer normal is ill-defined. To fix our ideas, we assume that x is an inner point of 
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the face 


7Ti := {a: 1 = 0} D dfl ; (4.26) 

let n(x) denote the unit outer normal at x, so that n(x) = (—1, 0,0). On the 
one hand, we can define 


P r n (x) := RP 


<0|7™(x)|0Mx) = -RP (0|^( X )|0) . 


(4.27) 


On the other hand, we can consider the alternative definition 


pr(x) 


lim (0|T ij (x / )|0) re n)^(x) 

X u.X ^X 


lim (0|Tji(x / )|0) r 

x'eO,x'-s>x 


(4.28) 


As a matter of fact, the alternatives 114.271) 114.28)) give the same result for the renor¬ 
malized pressure ; the rest of the present subsection is mainly devoted to the justifi¬ 
cation of this statement, which requires a nontrivial analysis. 

Consider the expression (I4.24p of (0|Tj“|0). Due to the considerations in the previous 
subsections, it is apparent that the terms in (j4.24j) deserving special attention when 
comparing the definitions (I4.27j) (I4.28p for the pressure at x = (x 1; x 2 ) G are 
those with n, cr such that 


^Lr( x i, xi) = 0 and 9t£ (xi, xj) ^ 0 for x' = (xj, x' 2 ) G Q 


these are easily seen to correspond to the choices 


(4.29) 


n = d\ — 1 and o G such that cr(di) = 1 . 


(4.30) 


Indeed, all the terms corresponding to values of n, u different from the above ones are 
straightforwardly seen to yield the same results according to both the prescriptions 
(14.27P and (14.281) . Now, let us focus on the potentially troublesome terms described 
in Eq. (14.241) . corresponding to a choice of the form (14.300 ; in the sequel we will 
show that these terms do not contribute to p[ en (x) for both the alternatives (14.2 71) 
f]4.28p . In order to prove this, let us denote with x' = (xj, xj) either the previously 
mentioned boundary point x G 7Ti or a point of fh In the expression (I4.24p for 
(OlT^x'^O) we pick up any problematic term with n, a as in Eq. (I4.30P : this reads 

+ + g^-d+i (m 2< yi 2 dl -i,a)') + ( 4 -31) 


+ m 2 f - ( \ Snd^dy't + (\~Z) d x n y n - f d x n x n J 0^ 


d+1 


(m 2 K- 1,.) 
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(some of the arguments have been suppressed for the sake of brevity). With some 
effort, noting that 9 r $ 1 _ 1)ff (x'i,y' 1 ) = (d[+y'\) 2 + Y^Wi-y*i ) 2 (duetoa(di) = 1) 
and using Eq. (14.151) , we can re-write expression (14.31 j> as 


5il 1 4 ^ 


(d+1 — u)& u-d-i ( z 2 ) + (5 u~d+ 1 ( z 2 ) — z 2 &u-d~ 


3 (Z 


z=2mx'\ ( 4 . 32 ) 


We now claim that 

f(u, 0) 


and 

lirn 

b'J-K) \ 

let us remark that in both the above equations the prescription of taking the regular 
part is superfluous, and | u=0 indicates the analytic continuation at u — (1. Eq.s (14.33j) 
(14.34)1 state that the problematic terms do not contribute to p[ en (x) as defined by 
Eq.s (I4.27P and (14.28(1 . respectively. 

In order to prove Eq. (14.3311 we note that 


= 0 


u =0 


= 0 ; 


(4.33) 

(4.34) 


f(u, 0) 



d+l—u 

2 


u - 


-d — 1 


+ r 


u - 


-d+1 


0 (4.35) 


where in the first passage we used Eq. (14.1611 . while equality to zero follows from 
the well-known relation r(^+l) = zT(z). 

To prove Eq. (14.3411 . recalling the definition (j4 .13 [1 we infer (for all x'\ > 0) 


f(u,x'\) 


IL =0 


(4.36) 




0 2 f(d+l)Kd+i (z) + zKd-i (z) — z Kd +3 (^)") 

V 2 2 2 / _ 


z=2mx'\ 


= 0 


in this case equality to zero follows from the identity below, holding for Bessel 
functions K u of any order (see |27], p.251, Eq.10.29.1): 


zK v+1 (z) - zK u _i(z) = 2v K v {z) . (4.37) 

In the above we assumed x to belong to the face with x l = 0 but, of course, similar 
considerations also hold for all the other boundary points not on the corners. 

Now, let us spend a few words about points on the corners of dQ, which appear if 
d\ > 1; we already noticed that the outer normal is ill-defined at these points, so 
that the notion of pressure is itself problematic. The natural strategies that could 


33 






































be guessed to overcome the problem make apparent some pathologies that we prefer 
to describe in an example, rather than in general: see subsection 14.91 

In passing, let us anticipate that an analysis similar to the one of this subsection will 
also be given in Part IV (see subsection 3.6 therein), for the case of a massless field 
confined within a d- dimensional box and fulfilling Dirichlet boundary conditions. As 
in the present setting, the alternative definitions (12. 5 7ft (I2.58P will be found to agree 
at all boundary points except those on the corners, where pathologies appear. 

4.7 Introducing two examples. The framework developed in the previous sub¬ 
sections will be illustrated heafter, for d = 3, in these cases: II := (0, +oo) x R 2 , 
representing a half-space, and hi := (0,+oo) 2 x R, representing a wedge bounded 
by orthogonal half-planes. In both cases, we consider Dirichlet and/or Neumann 
boundary conditions. 

4.8 A half-space in spatial dimension d — 3. Let 

Q := (0, +oo) x R 2 ; (4.38) 

this is the subcase of the general setting (14. ip corresponding to d = 3 and 

dx = 1 , d 2 = 2 . (4.39) 


With the above choices, the symmetric group appearing in the general framework 
of subsection 14.51 consists of the sole idendity ( S ^ — Si — {id}). We have X] = 
(a^) = x 1 (and the analogous relations for y x ); besides, OTo^xi, y x ) = \x 1 +y l \ and 
W Md (xi,yi) = | x 1 —y 1 \ . Using the relations (l4.17H4.19p . (I2.16H2.18P and Eq.s (I4.15P 
(I4.16P and (14.371) . with some simple algebraic manipulations we obtain the following 
expressions for the non-vanishing components of the regularized stress-energy VEV 
(where x = (x 1 , x 2 , x 3 )): 


(orewio) = - 


m 


rm\~ u 


32vr 3 / 2 r(H±±) \ K ) 


(l—w) r 


u — 4 


+ 


- 2 5 -! Ql ( ( i - n<2W^) + ( ^ - (3-u)a<5^(z 2 ) 


z=2mx 1 


<0|3Ti(x)|0) = 


<0|IS(x)|0> = (0|T 3 « 3 (x)|0) = 


mn 


32^3/2 r(H±i) \ K J 


m 


32v r 3 /2r(H±l) V K 


2 5 "2 Ql ((:-n(g^) + 


2 

2 — u 


( n> ^ u v ( 11 —^ 
u —4 


m 


+ 


(3-u)^©^(2: 2 ) 


z=2mx 1 . 


(4.40) 

(4.41) 

(4.42) 
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The above expressions are easily seen to give the meromorphic continuation of the 
regularized stress-energy VEV to the whole complex plane, with poles determined 
by terms with the gamma function. In particular, all the above components have a 
simple pole in u — 0; thus, we follow the extended version of the zeta approach and 
define the renormalized quantities to be the regular parts in u — 0. Recalling again 
that Eq. (I2.14p gives 



we write the final results in the form (I2.13p . obtaining 

(0|ho(x)|OU^(3(41 n (B) +1 ) +32Ql «) + 

/ 1\ m 4 / (2mx 1 )K 1 (2mx l ) + 3K 2 (2mx 1 ) 

1 y 6 / 7T 2 y (2 mx 1 ) 2 

<0|?„(x)|0U = (41„(^) - 3) ; (4.44) 

<0|f 22 (x) |0}„„ = (0|f s3 (x) |0>_= - ( 3 (41„ 0 - 3) + 32 a, + 

/ 1\ m 4 f {2mx 1 )K l {2mx 1 ) + 3K 2 (2mx 1 ) 

1 \ 6/ 7T 2 y (2'ma: 1 ) 2 

Let us comment briefly on the above results. Firtly, note that the renormalized VEV 
of the stress-energy tensor does not depend explicitly on the spatial coordinates 
x 2 ,x 3 ; this comes as no surprise due to the homogeneity with respect to these 
variables of the spatial configuration considered. Besides, in agreement with the 
general results of subsection 14.61 the components (0|T ] ! j|0), (0|Tii|0) ren are constant 
and the two alternative definitions for the pressure on the plane Hi = {x\ = 0} (see 
Eq.s (14.2711 (I4.28P ) give the same result; more esplicitly, we obtain 

vT = —M 0\T u \0) ren (1 e {1,2,3}) . (4.46) 

In conclusion, let us make a comparison with the results derived in [2] for the 
configuration with a single plane in arbitrary spatial dimension (to be considered 
here with d — 3); therein the attention is restricted to the “minimal” (£ = 0) and 
“conformal” (£ = 1/6) settings. In both cases the results derived here are found to 
agree with those reported in [2] (let us remark that the mass scale k employed here 
does not coincide with the one considered therein; the latter is proportional, via a 
numerical coefficient, to m 2 / k). 


(4.45) 
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4.9 The rectangular wedge. Let us pass to the case of a wedge in R 3 , bounded 
by two perpendicular half-planes; this is represented as 

n := (0, +oo) 2 x R , (4.47) 

corresponding to the general framework (14.11) with d — 3 and 

d 1 = 2 , d 2 = 1 . (4.48) 

In passing, let us mention that the rectangular wedge model is also considered by 
Actor and Bender |2], for arbitrary spatial dimension; yet, these author restrict the 
attention to a massless field, a case we discuss in the next subsection. 

In the present setting, the symmetric group (S^ = S 2 ) of subsection 14.51 consists of 
two elements, i.e., the identity id and the exchange p: 


S dl =S 2 = {id,p} , id( 1) = 1, id( 2) =2, p(l) = 2 , p(2) = 1 . (4.49) 

Moreover, with Xl = {x\,xf) = (ad, a; 2 ) and the analogous relations for y 1; we have 

9V( Xl , yi ) = ((x 1 +y 1 ) 2 + (x 2 +y 2 ) 2 V 

) (i /2 for a e S 2 = {id,p} ; (4.50) 

9Wx 1>y i)=({x 1 -y 1 f+{x 2 -y 2 f 


/ \l/2 / \1/2 

*Mx 1>yi ) = ((z 1 -? 1 ) 2 ^- {x 2 +y 2 f) , = ((ad+y 1 ) 2 ^- ( x 2 -y 2 ) 2 ) . 

Also in this case, we can use the relations (12.10112.151) . (j4.17H4.lHp and the identities 
(14.15P (14.16p and (I4.37P to deduce expressions for the non-vanishing components of 
the regularized stress-energy VEV. More precisely, we obtain the following (with 
x = (ad, x 2 , x 3 )) 


<0fe(x)|0) = - 

- 2 5_i E ° 


m 


m 


(i—«) r( 


u—4 


r U 
2 5 "2 


327T 3 / 2 r(^±i) V K 

(3-u)^j&u^{z 2 ) + - ^J<8h=2(z 2 ) 

« ia 2 - (2-u)^0u-4(z 2 ) +Q ^^®u=2(z 2 ) 


+ 


(4.51) 


2 = 1.2 


+ 


z=2mx x 


z=2m^J (x 1 ) 2 +(x 2 ) 2 


36 



















(0|T”(x)|0) = <0|T“(x)|0) = 


m 


32 vr 3 /2r(“±i) \ K ) 


m \-u 


r( + (4-52) 


- 2 5 * a pU) Stj ~ ( 3 -u)^)&u^(z 2 ) + (^-^)&u^{z 2 ) 

2 5 ~% a x a 2 (<% ^(3-M)0iizd(^ 2 ) + (3u^2 (z 2 )j + 


+ 


z=2mxVU) 


4 m 2 x' l x j 6.-6 (z 2 , . , _ 

2 Jz=2myJ ( x 1 ) 2 +(x 2 ) 2 


(0|r 3 " 3 (x)|0) = 


m 


fmy 

32 ^ 3 / 2 r(H±i) 


■u —4 


+ 


for i,j G {1,2}; 

(4.53) 


- 2 5 -* E 


a,; 


»= 1,2 


- (3-w)^© 2f 4(^ 2 ) + Q - ^)®^2{z 2 ) 


+ 


z=2mx l 


- 2 5 2 ai a 2 (( - (2-M)()®-i(z 2 ) + (^- ^j&u-2(z 2 


' z=2myj (x 1 ) 2 + (x 2 ) 2 

The above expressions give the meromorphic continuation of the regularized stress- 
energy VEV to the whole complex plane, with poles determined by terms with the 
gamma function. In particular all components (14. 5 1114.531) have a simple pole in 
u — 0, and we must resort to the extended zeta approach taking again the regular 
parts at u — 0. Recalling once more that Eq. (\2.14\ gives 


£3 r 1 

0 

hereafter we report separately the conformal and non-conformal parts of each com¬ 
ponent of the renormalized stress-energy VEV; these are 


(0| T oo(x)|0) re n = 


mn 


384vr 2 


3(4, n (£) +1 ) + 32 E 


.2k, 

32 ol\Ol 2 


Otj 


i= 1,2 

zK^-K^z) 


Ki(z) 


+ 


z=2mx l 


m 


(01-^00( X ) |0)ren — o 


7T- 


ai 

.*= 1,2 


zR\(z) + 3 K 2 (z) 


z=2m v /(i 1 ) 2 +(i 2 ) 2 


+ 


(4.54) 


+ Oi\Oi 2 


z=2mx % 

z K 1 (z) + 2K 2 (z) \ 

z2 Jz=2m y /(x 1 ) 2 +(x 2 ) 2 


(4.55) 
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(0|f’»(x)|0>„„ = 


32 ol\OL2 


( 0 |^(x)| 0 ) ren = 


m 


384vr 2 


m 


3 ( 4 ln (—) — 3 ) + 32 a p(j) S %J 


Ki(z) 


+ 


Z Jz=2mxVU) 


Z J z=2myj (x 1 ) 2 +(x 2 ) 2 

for i,j G {1,2} ; 
zK 1 (z)+3K 2 (z) 


m 


7T 


a p(j) ^k 


+ 


z=2mxP(i) 


( s z Ki(z)+3K 2 (z) a 2 4 , K 3 (z) 

+ «ia 2 Oij ---4 m x x J — 5 — 

\ z 2 z s 

for i,j G {1,2} ; 


z=2myj (a; 1 ) 2 +(a; 2 ) 2 


(4.56) 


(4.57) 


I-^33( X ) 10)ren — 


mn 


384vr 2 


H 41n U- 3 ) +32 5> 

i=l,2 

zK!(z)-K 2 (z) 


M z ) 


+ 


z=2mx l 


32 CI 1 CI 2 


z=2myj {x 1 ) 2 -\-{x 2 ) 2 


(4.58) 


m 


I -^33( X ) 10)ren — n 


n* 


ai 

_ 2 = 1,2 


z K\(z) + 3/i2(^) 


+ 


+ 0i\0i 2 


z=2mx z 

zK 1 (z)+2K 2 (z) 


(4.59) 


Z J z=2m- s J (x^Y+ix 2 ) 2 _ 

As was to be expected, since the configuration (14. 4 71) is invariant under translation 
along the x 3 direction, none of the expressions (14.54114.59} depends on the spatial 
coordinate x 3 . 


Let us now discuss the pressure at points in the half-plane tti = {x 1 = 0, x 2 > 0}; 
note that, we are excluding points on the axis ( := {x 1= x 2 = 0}. The two alternative 
definitions (14.27} (14.28} are easily seen to give the same result for the renormalized 
version of this quantity, in agreement with the general results of subsection 14.61 
indeed, we can equivalently put x 1 = 0 in Eq. (14.52} for (0|T)“(x)|0) (i G {1,2,3}) 
and then analytically continue up to u — 0 , or directly evaluate the renormalized 
expressions 04.56} 04.57} for (0|Tb(x)|0) ren in x 1 = 0. I 11 both ways, we obtain 




— &il 


m 


384tt 2 


3 4 In 



—— (l + «i)a2 


( S )- 3 ) +32 “ 2 

z K 1 (z) + 3K 2 (z) 


{l + ai-L)zKi{z) + 3aiK2{z)\ 

i 2 ) + 

(iG {1,2,3}) . (4.60) 

z=2mx 2 
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Let us stress that, differently from all the configurations considered so far within this 
paper, in this case the renormalized pressure on the boundary depends in general 
on the parameter £; more precisely, this happens whenever Neumann boundary 
conditions are imposed on the half-plane 7iq (so that Oi = +1). 

Now, let us consider the axis ( = {a: 1 =x 2 = 0}; at any point of this axis the outer 
normal is ill defined, so that there is a basic obstruction to speaking of the pressure, 
ffowever, we can discuss what happens if a point x = (0, x 2 , x 3 ) G 7Ti moves towards 
the axis (, i.e., if we consider the limit x 2 —> 0 + . In this limit p[ en (x) is found 
to diverge; more precisely, Eq. (14.60P and the known asymptotic behaviour of the 
Bessel function K u near zero (see [2?], p.252, Eq.10.30.2) give the following: 


pT(x) = o 



for x G 7Ti, x 2 —y 0 + . 


(4.61) 


4.10 The previous examples in the zero mass limit. Let us first consider 
the half-space configuration (14.38(1 bounded by the plane 7r 1; analysed in subsection 
14.81 The expressions (I4.43H4.451) for the components of (0|T Atl/ |0) ren give, in the limit 
m —* 0 + , 


(0|T^(x)|0) ren 


/-I 0 0 0\ 

/ 1\ 3 aq 0 0 0 0 

\ 6/ 87r 2 (x 1 ) 4 0 0 1 0 

\ 0 0 0 1 / 


(4.62) 


as for the pressure on a point x G 7Ti, starting with Eq. (|4.46[) and taking the limit 
m —y 0 + , it is trivial to infer 


r(x) = 0 (i G {1,2,3}). (4.63) 

Let us pass to the case of the rectangular wedge (see subsection 14.91) . Eq.s (14.541 
14.591) for the renormalized stress-energy VEV give, in the limit m —y 0 + , 


(0|f M „(x)|0) ren 


/x,i/=0,1,2,3 


/-1 0 0 o\ 

oio 2 0 4i(x) -B(x) 0 

967r 2 p 4 0 f?(x) 4 2 (x) 0 

\ 0 0 0 1 / 


(4.64) 




( 

-1 

0 

0 

°\ 


( 

-C'o(x) 0 0 

0 ^ 


(-3 

CM.lCy.2 


0 

4i(x) 

B(x) 

0 

3 


0 Ci(x) 0 

0 


87i 2 p 4 


0 

B(x) 

A (x) 

0 

8n 2 


0 0 C 2 (x) 

0 




V 

0 

0 

0 

1 ) 


V 

0 0 0 

Co ( x ) j 



5(x) 


A(x) := 
Ax'-x 2 



Q(x) 


a p(0 

(x p b )) 4 


C 0 (x) := 


OL 1 


[X 


1)4 


+ 


Qy 


[x 


2U 


0:10:2 


for i — 1, 2 ; 
p := a/( x 1 ) 2 + (a: 2 ) 2 
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(recall that p(l) = 2, p(2) = 1). Next, consider the expression (I4.60p for the 
pressure acting on a point x in the half-plane i\\ = {x l = 0, x 2 > 0}, for m > 0; 
using the asymptotic behaviour of the Bessel function K v near zero (see \27\, p.252, 
Eq.10.30.2) we infer, in the limit m —>■ 0 + , 



(i e {1,2,3}) . (4.65) 


Notice that, as for the massive analogue (I4.60jh the above expression for the renor¬ 
malized pressure depends explicitly on £ if we assume Neumann boundary conditions 
on 7Ti (so that aq = +1). 

In passing, let us remark that both results (I4.63P and (I4.65P for the pressure on 
7T\ could be determined equivalently via the prescription (I2.58P : according to the 
latter, we should have first considered the renormalized stress-energy VEV inside 
the corresponding spatial domain in the limit of zero mass (see Eq.s (I4.62p (14.641) 1 
and then move to the boundary (half-)plane 7Ti, i.e., take the limit x 1 —>■ 0 + (3)- 

Let us comment briefly on the construction described above, namely, that of taking 
the zero mass limit (m —>■ 0 + ) of the renormalized results (I4.43M.45P (I4.46P and 
(I4.54M.59P (14.601) for the theory with a massive field. As a matter of fact, this 
procedure corresponds to studying the case of a massless scalar field (in the same 
spatial configurations) with the technique of subsection 12. 141 Indeed, in the massless 
case (m = 0) the spectrum of the fundamental operator A = —A is [0, +oo), for both 
the settings (I4.38P and (14.471) : according to the framework of the cited subsection, we 
could treat these cases using the deformed operator A e := A + £ 2 (see Eq. (5.14)), 
and eventually taking the limit £ —> 0 + . On the other hand, if e is identified with 
m, we recover the present constructions (pt). 

Summing up: Eq.s (I4.62M.65P and (I4.64M.65P yield the renormalized VEVs of the 
stress-energy tensor and pressure for a d = 3 massless held, respectively confined 
within a half-space and a rectangular wedge, fulfilling either Dirichlet (cq = —1) or 
Neumann (cq = +1) boundary conditions. 

We point out that, in the Dirichlet case, the above results are found to agree with 
those derived by Actor and Bender [2] for £ = 0 and £ = 1/6, via a different 
version of the zeta approach (also involving, essentially, a subtraction of divergent 
contributions; see subsections 3.1.1 and 4.1 of the cited paper, setting d — 3 therein). 

12 Notice that not all of the components in Eq.s (14.621) (14.641) are finite for x 1 —> 0 + , but only 
those involved in the computation of the pressure on 7Ti ; for example, in both cases we have 


a; 1 —»0+ 


lim (0|f 22 (x)|0) 


ren 


= OO . 


13 The same comments could be made in general for any spatial dimension d and for any number 
d\ of faces. 
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Let us also mention that the massless half-space configuration in the case of Dirichlet 
boundary conditions was considered as well in our previous work HZ|; therein the 
same results were obtained starting from the renormalized stress-energy VEV of 
a massless held between parallel planes, and taking the limit of infinite distance 
between the planes (see Eq. (5.7) in [IT]; after an exchange 1 f> 3 in the coordinate 
labels this becomes the present Eq. (14.6211 . with aq = —1). 

In the next section we show that the renormalized expressions (14.62114.67)1) . here 
deduced as the zero mass limit of a massive theory, can be obtained equivalently as 
particular cases of the theory of a massless scalar held confined within two half-planes 
forming an angle a of arbitrary width; more precisely, the present configurations with 
one single plane and two orthogonal planes correspond, respectively, to the limits 
a —> Ti and a — > 7t/2 . 

5 The case of massless field in a 3-dimensional 
wedge 

5.1 Introducing the problem for arbitrary boundary conditions. In this 
section we consider the case of a scalar held (with no external forces) confined 
within a 3-dimensional wedge, meaning that the spatial domain is the portion of 
the space R 3 enclosed by two half-planes 7r 0 , 7T a forming an angle a 6 (0, 27r] ; we 
assume 7T 0 = {x 2 = 0, x 1 ^0}. Suitable boundary conditions will be specihed in the 
following. 

We choose to confine our attention to the massless case (V = 0) since, in this 
case, we are able to perform the explicit computations with a moderate effort for 
arbitrary values of the angle a e (0, 27r]; in the massive case (V = m 2 ) we could 
give an exhaustive analysis for rational values of a/n but this would require a big 
computational effort and produce cumbersome expressions for the final results. 

We already pointed out in the Introduction that the present framework has been 
previously considered by Dowker H21Q3], Deutsch and Candelas HU, Saharian et 
al. [2H) SIJ and by Fulling et al. [23]) amongst other; a more detailed comparison 
between our results and these works will be performed in subsections 15.6115.91 
In passing let us notice that, in case of either Dirichlet or Neumann boundary 
conditions, the spatial domain under analysis corresponds for a = n and a — 7r/2, 
respectively, to the configurations with a boundary made of a single plane and of two 
orthogonal half-planes. Besides, let us also stress that for a = 2i r the two half-planes 
7To, n a overlap; because of this, for Dirichlet or Neumann boundary conditions, the 
boundary consists of a single half-plane, while, in the case of periodic boundary 
conditions, the spatial domain can be identified with R 3 . So, in the latter case 
one is actually considering a massless scalar held on the whole Minkowski spacetime. 
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We will comment further on each specific case in the next subsections. In particular, 
we will show that in the cases with a = n and a = n/2 one recovers the results of 
subsection 14.101 

In order to deal with the present configuration, it is advisable to pass to the system 
of cylindrical coordinates 

x (->■ q(x) = (p(x), 6*(x), z(x)) G (0, +oo) x (0, 27r) x R , (5.1) 

whose inverse will be written q i—>• x(q). These coordinates are chosen so that 
z(x) = x 3 and the boundary dfl corresponds to the limit values 6 = 0 and 6 = a\ 
the spatial line element reads 

cU 2 = dp 2 + p 2 d6 2 + dz 2 . (5.2) 

In order to avoid clumsy notations, given any function —> Y, x H > /(x) (with Y 

any set), the composition q G q(h2) /(x(q)) will be written as q H> /(q). 

Since curvilinear coordinates are being employed, one should refer to the framework 
of subsection 12.151 Writing p, 9 ,2 for the coordinate labels, we see that the only 
non-vanishing Christoffel symbols associated to the line element (15.2ft are 7 = —p 
and 7 e pd = 7 e dp = i; so, the second order covariant derivatives of any scalar function 
/ are given by 

^ppf = d pp f , V p ef = dpgf - i d e f , V pz f = d pz f , 

Veef = deef + pd p f , V 0z f = d ()z f , V zz f = d zz f . 

In conclusion, let us stress that the configuration under analysis could be dealt with 
as a slab configuration where = fix x R, and Qi C R 2 corresponds to (0, +cxd) x 
( 0, a) in terms of the coordinates ( p,6 ); yet, this approach is not convenient. In 
fact, if one works directly on the 3-dimensional spatial domain hi (for any of the 
boundary conditions to be considered in the following), the modified cylinder kernel 
T defined in Eq. (I2.28[) associated to the fundamental operator A = —A can be 
expressed in terms of elementary functions (moreover, it is meromorphic in the 
variable t on the whole complex plane). On the contrary, to treat the problem as 
a slab configuration we should use the analogous kernel T for the reduced operator 
A\ on Oi, or any other integral kernel related to the latter; these kernels do not 
possess simple expressions (only integral representations are available), so that the 
whole analysis would become a lot more involved. 

5.2 The Dirichlet kernel. For any of the boundary conditions to be considered 
in the following, it turns out that the spectrum of A is [ 0 , + 00 ). Since { 0 } is a 
non-isolated point of the spectrum, we must resort to the methods discussed in sub¬ 
section 12.141 to determine the renormalized Dirichlet kernel (and its derivatives); we 
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regularize the theory using the deformed operator A e = ( \fA-\-e) 2 (see Eq. (5.15)), 
whose choice is found, a posteriori, to be more effective from the computational 
viewpoint. In the sequel (see subsections 15.6115.9p we derive the explicit expression 
of the modified cylinder kernel T of A, for several types of boundary conditions; in 
any case T is found to be a meromorphic function of t, decreasing faster than tr 1 for 
9?t —* +oo; this allows us to proceed as explained in subsection 12. 141 In this way we 
obtain, for the Dirichlet kernels and its derivatives, these renormalized expressions 
at s = ±1/2, respectively (see Eq. (I2.73P ): 


£>i(q,p) = Res(2t 3 T(t;q,p; 


p=q 


;0 


p—q 


p—q 


V 2lu £) + i (q, p) = Resit ^^(t; q, p 


;0 


p=q 


for v,w any two cylindrical coordinates . 

Here and in the sequel, we are using the following notations: 

q=(p,e,z), |> = ■ 

/(x(q),x(p)) = /(q,p) for any /:Cx!f-l Y (Y a set) 


(5.4) 

(5.5) 


(5.6) 


5.3 The stress-energy tensor. Relations (I5.4j) (15.5p . along with Eq.s (l2.16H2.lgp . 
can be used to obtain the following expressions for the renormalized VEV of the 
stress-energy tensor: 


(0|Too(q) |0)ren = ReS 


(‘" 3 [(: 


+ 2£) T(t; q, p) + 


± () t d q dpeT( t;q,p) 


;0 ; 


J p=q 


(5.7) 


(0|±ij(q)|0)ren = (0|T j? ;(q)|0) ren = Res ( t 3 ( ^ - 2^) % T(t;q,p) + 


t 2 ( ( \ ~ \^ 1 -£) Sijd^dpi ) T(t; q, p) 


;0 


J p=q 


for i,je{p,e,z] . 


(5.8) 


In the following subsections we will first compute T and then explicitly evaluate the 
residues appearing in Eq.s (15.7p (15.8p for the cases where either Dirichlet, Neumann 
or periodic boundary conditions are prescribed; we will present the final results in 
the form (12.1311 . recalling once again that Eq. (12. 14p implies 


£s = 


1 

6 ' 


(5.9) 
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5.4 The boundary forces. In order to discuss this topic, since the spectrum of 
A contains {0} as a non-isolated point, we must resort once more to the framework 
of subsection 12 . 141 again, we choose to use the deformed operator A e := (v^A + s) 2 - 
We can consider the two alternative definitions (I2.66P (12. 6 7ft for the renormalized 
pressure on the boundary. Contrary to the configurations analized in the previous 
sections, in this case the two mentioned alternatives do not yield the same result. 
As an example, let us focus on the pressure acting on the half-plane 7 t q . We indi¬ 
cate with n(q) the unit outer normal at a point of this half-plane with coordinates 
q = (p, cx,z); this has components (n p ,n e ,n z )( q) = (0,1/p, 0). On the one hand, 
prescription (I2.66p corresponds to put, for i G {p, 9, z} , 


P? n { q) : 


lim 

£- 5 - 0 + 



0 (0|T’“(q)|OK(q)) 


1 

P 


lim RP 

£—> 0 + 


O (0re“(q)|0) 


(5.10) 


(first consider the regularized stress-energy VEV on 7r Q , then analytically continue 
at u = 0 and finally take the limit e —>• 0 + ). Similarly to Eq. (15. 8[ . the prescription 
(15T0D yields 


p? n ( q) 


- Res 3 



2f) S i0 T( t;q,p) + 





T( t;q.p) 


;0 


p=qS7r a 


(5.11) 


let us stress that in the above expression, we first perform the evaluation on the 
boundary and then compute the residue. 

On the other hand, the alternative prescription (12.67ft yields, for i G {p, 9, z } , 


P? n ( q) : = ( , lim (0|Tij(q / )|0)ren n J (q) = 

.q'vq,x(q')ef! 


. , lim (0|Ti. e (q')|0) r 

P Vq , ^-q.x(q')eO 


(5.12) 


In the cases to be considered in subsections 15.6115.81 it will be apparent that the ex¬ 
plicit expressions obtained for the renormalized stress-energy VEV inside the wedge 
diverge when approaching the boundary, in such a way to make divergent the renor¬ 
malized pressure defined by 115. 12)1 . Because of this, we will always refer to Eq. 
(15.lip to deal with the pressure on the boundary. 

5.5 Some remarks. For the computation of T we will often refer to [23], where 
this kernel was already determined for the present configuration (but used in a 
different renormalization scheme, based on point splitting as an alternative to zeta 
regularization; we note that the kernel T in [23] is the opposite of the kernel T 
considered here). 
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At the end of each subsection dealing with Dirichlet and/or Neumann boundary 
conditions, we will comment briefly on the results obtained for the renormalized 
stress-energy VEV (l5.7H5.8p and pressure (I5.12p in the cases a = tt and a = ir/2, 
respectively, describing a massless scalar field on a half-space and inside a wedge with 
orthogonal half-planes. Recall that these very same configurations were analysed as 
the zero mass limit of a corresponding massive theory in subsection 14.101 indeed, we 
will End that the same results obtained therein can be re-obtained from Eq.s (15.71 
15.8p and (I5.12p . returning to the Cartesian coordinates x 1 =_psind, x 2 = p cos 9, 
x 3 = z and considering the appropriate transformation laws (l 14 l ). 

5.6 Dirichlet boundary conditions. Let us first consider the case where the 
held fulfills Dirichlet boundary conditions on the half-planes 7To, 7i a . In this case 
a complete orthonormal system of (improper) eigenfunctions (F k ) k £K of A = —A, 
with eigenvalues (a J k )k&K.-, is given by 

-Ffc(q) := 

l o\ := a; 2 + h 2 for k = (n,u, h) e /C = N x (0, +oo) xR 


'—J Xn (ujp) sin(A n 9)e 
na 


ihz 


TiTT 

An • j 

a 


(5.13) 


(here and elsewhere we are considering the set of positive integers with N : = 
{1,2,3,...}; /C is endowed with the counting measure on N times the standard 
Lebesgue measure on (0,+oo) x R, meaning that J K dk = Xln=L ./ 0 +O ° du dh). 
The modified cylinder kernel T can be evaluated starting from its eigenfunction 
expansion (12.281 . which in the present setting reads 

^ r p tu k _ 

T( t;q,p)= / dk - F k (q)F k (p) = (5.14) 

JK u k 

i +°° r+oo r+oo e ~Wu‘ 2 +h 2 

— y>n(Anfl) sin(Anfl') / du: u J Xn (cop) J\ n (cup') / dh—=== e lh{z ~ z ' ] . 

7o J -oo Vco 2 + h 2 

With some effort, the integrals over h and oj in the above expression can be explicitly 
evaluated to yield 


+oo 

f(t;q,p) =- r —— V]sin(A„6>) sin(A„6> / ) e~ XnV , 

7rex pp smh v 

n=1 (5.15) 

v := — In (— -— J , r± := y /Tp± p ') 2 + (z — z') 2 + t 2 . 

\r + + r_ J 

We refer to [22j for the detailed computations giving the above result; the notations 
v, r± are mutuated from this reference (also see [23]). The series in Eq. (15. 15[) can be 

14 Clearly, (0|T , Ml/ (q)|0) T . en and p^ en { q) transform, respectively, as a rank-two tensor and a vector. 
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re-expressed via four geometric series writing the trigonometric functions in terms 
of complex exponentials; in this way we obtain 


T( t;q,p) = 

1 / cos (^(9 — 6')) — e~^ v cos (^(6 + 6')) — e~^ v \ (5.16) 

47TQ! pp' sinh v 1 cosh(^n) — cos (^{9 — 6')) cosh(^u) — cos(^(6*+6 1 ')) J 

Now, we resort to Eq.s (j5.7H5.8p for the renormalized stress-energy VEV; evaluating 
the residues therein we obtain 


(0|f^(q)|0) re n 


A(q) 


f-1 0 
0 1 
0 0 
\ 0 0 


0 

°\ 


/-(£(q) + C(q)) 0 0 0 \ 

0 

0 

-H) 

0 5(q) -pE( q) 0 

3p 2 

0 

0 -pE( q) p 2 C( q) 0 

0 

V 


\ 0 0 0 5(q) + C(q)/ 


A(q):= 


7r 


-a 


B(q) := 


97T 4 —3^ 2 (27r 2 +a 2 ) sin 2 (—)+a 2 (;r 2 —a 2 ) sin 4 (—) 


14407 r 2 a; 4 p 4 ’ 

/ 37T 2 — (7T 2 — a 2 ) sin 4 (—) 

C{ q) := . 2 \ . 2 L. \ aJ 

87r 2 a 2 sm (—) p 4 


247r 2 a 4 sin 4 (—) p 4 


E( q) := 


3t r cos(^) 
8a 3 sin 3 (—) p 


(5.17) 


It can be easily checked that the above result agrees with the one derived by Saharian 
and Tarloyan by means of point-splitting regularization in [31] (see, in particular, 
Section 3 therein); see also the former papers by Dowker et al. |12, [13] and by 
Deutsch and Candelas HU, where point-splitting regularization is used for the com¬ 
putation of the conformal stress-energy VEV alone. 

Let us briefly comment on the explicit expression (15.171) obtained for the renor¬ 
malized VEV (0|T /w (q)|0) ren . First of all, notice that the function A(q) multipling 
the conformal part of the renormalized VEV is positive for a < 7r, negative for 
a > 7i and vanishes for a = n. Next, let us stress that both the conformal and 
non-conformal parts diverge quartically in p in the proximity of the axis {p = 0}. 
The non-conformal part also diverges near the half-planes 7To, 7 t q , that is for 9 —y 0, a; 
because of this, the pressure on these half-planes evaluated according to Eq. (15. 12)1 
is infinite. On the other hand, the alternative definition (I5.10H5.1T]) (first move to 
the boundary, and then take the analytic continuation) gives a finite pressure on n a 
with components 


P? n ( q) = -$ 


4 4 

7T — a 


iO 


4807r 2 a 4 p 3 


(5.18) 


To conclude, we consider the special cases a = tt and a = n/2 (a space domain 
bounded by a plane or by two perpendicular half-planes), and compare the present 
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results with the ones of subsection 14.101 This can be done with the procedure out¬ 
lined in subsection 15.51 (i.e.. returning to Cartesian coordinates via the appropriate 
transformation rules for tensor coefficients). Indeed, the expressions (15.171) (15. 18[) 
(giving the renormalized stress-energy VEV and pressure) are easily seen to yield 
for a = 7T and a = tt/2, respectively, Eq.s (I4.62|) (j4.63j) with aq = —1, and Eq.s 
(I4.64p (I4.65P with aq = a 2 = —1. 

5.7 Dirichlet-Neumann boundary conditions. Let us now pass to the analysis 
of the wedge configuration where Dirichlet and Neumann boundary conditions are 
prescribed, respectively, on the half-planes 7r 0 and n a ; let us point out that, to the 
best of our knowledge, this case has never been considered before in the literature. 
A complete orthonormal system of (improper) eigenfunctions (Ea,)*,^ of the funda¬ 
mental operator A, with eigenvalues (ca 1 2 )^*;, is given by 

Efe(q) : = 

u\ := oj 2 + h 2 for k = (n, c u, h) G /C = N 0 x (0, +oo) x R 



J Xn (up) sm(X n e) e 


ihz 


\ n := n + 


7T 

a 


(5.19) 


(N 0 := {0,1, 2,...} is the set of non-negative integers; the measure on the label space 
fC is such that J K dk = Yln=o / 0 + °° du dti). Resorting again to Eq. (12.28P and 
proceeding similarly to the case with Dirichlet boundary conditions, we can express 
the modified cylinder kernel as 

E(t;q,p) = 

1 / sinh(^u) cos(^-flQ) _ sinh(£tQcos(£(fl + flO) \ (5.20) 

2na pp/ sinh v \cosh(^u) — cos(^(0 — 6')) cosh(^u) — cos^^-t-fl')) / ’ 

where v si dehned as in Eq. (15.151) . Using Eq.s (15.7H5.8p we obtain the renormalized 

VEV of the stress-energy tensor: 


A(q) 


/I 0 

0 -1 


0 

\0 


0 

0 


0 0 
0 0 
3 p 2 0 
0 -1 


(0|TV(q)|0) r 


[i,i'=0,p,0,z 


/ 

1 \ 

-(?- 

~e) 


\ 

A( q) := 


0 
0 
0 

77r 4 + 8a 4 
11520vr 2 a 4 p‘ 


B( q) 

~PE( q) 
0 


,4 ’ 


0 0 \ 

~PE( q) 0 

p 2 C(q) 0 

0 B(q) + C(q)J 


B( q) : = 


-3vr 2 cos(^)(ll7T 2 —2a 2 + (vr 2 + 2a 2 )cos(^)) + 2a 2 (vr 2 + 2a 2 ) sin 4 (^) 

967r 2 a 4 sin 4 (—) p 4 


, , 67r 2 cos(—) + (7r 2 + 2a: 2 ) sin 2 (— ) / 3tt(3 + cos(^)) 

c <0 :=- a . ^(q) = oo • ( 5 - 21 ) 


167t 2 q; 2 sin 2 (—) p 4 


32« 3 sin 3 (—) p 4 
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Let us compare the above results with the ones of Eq. (I5.17p . holding for the case 
of Dirichlet boundary conditions on both the half-planes 7To,7r Q . As in Eq. (I5.17p . 
both the conformal and non-conformal parts of the renormalized stress-energy VEV 
diverge for p —> 0; the latter also diverges for 6 —> 0, a, so that Eq. (I5.12p yields 
again a divergent pressure on the boundary. On the other hand, the present results 
differ from the ones derived in the previous subsection because of some crucial 
features; in particular, the conformal part has an overall minus sign and the function 
A(q) in Eq. f!5.2ip is always strictly positive (whereas the one in Eq. (I5.17P changes 
sign for a < n and a > 7r). 

As for the boundary forces on 7T a , resorting to Eq. (15. lip , in this case we obtain 


77r 4 + 8a 4 
8tt 2 p 3 480a 4 


H) 


(5.22) 


We notice that also in this case the parameter £ appears in the final expression 
for the renormalized pressure; because of this the resulting boundary forces can be 
either attractive or repulsive, depending on the value of £. 

Again, we conclude comparing the results obtained for the renormalized stress- 
energy VEV and pressure for a = tt/2 with the analogous ones deduced in subsection 
14.101 Also this time, Eq.s (I5.2ip (15.221) (with a = 7t/ 2) are found to give, respectively, 
Eq.s (14.641) (14.65P with aq = 1, a 2 = — 1. 

5.8 Neumann boundary conditions. We now analyze the case where the field 
fulhlls Neumann boundary conditions on both the half-planes tiq, 7i a . A complete or¬ 
thonormal system of (improper) eigenfunctions {F k ) k ^ in L 2 (fl) of the fundamental 
operator A, with related eigenvalues (u>l) k& jc, is 





nn 


Fk( q) = \/—J\ n (ujp) cos(A n B) e ihz , A n := — , 

7ra a (5.23) 


col = lo 2 + hr 


for k = (■ n , u>, h) e /C = N 0 x (0, +cxd) x R 


(recall that N 0 := {0,1, 2,...}; again, we assume the measure on the label space /C is 
such that f K dk = /o + °° f- ^ dh)- ^-l so i n fhi s case, the modified cylinder 

kernel T can be evaluated according to Eq. (12.281) . More precisely, proceeding as we 
did in subsection l5.6l for the case of Dirichlet boundary conditions (see, in particular, 
the derivation of Eq. (15.16P ). we obtain 

t;q,p)= (5-24) 

1 ( - cos(|(g-gQ) e~ v - cos(^(g+gQ) 

Fttol pp' sinh v l cosh(^u) — cos(^(6 l — 9')) cosh(^u) — cos(^(6*+6 )/ )) 
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again, v si defined as in Eq. (I5.15I) . Now, we can resort once more to Eq.s (I5.7ti5.8p 
to evaluate the renormalized VEV of the stress-energy tensor; the result is 


(0|f #H ,(q)|0) ren 


A(q) 


/-I 
0 
0 
V 0 


0 0 
1 0 
0 —3 p 2 
0 0 


0\ 

0 

0 

V 


(5.25) 


+ 



/-(B(q) + C(q)) 0 0 0 \ 

0 B(q) —pE( q) 0 

0 -pE( q) p 2 C{ q) 0 

\ o 0 0 B(q)+C(q)J 


+ ^(q) 


/—2 0 0 0 \ 
0-100 
0 0 3p 2 0 

\ 0 0 0 2 


where the functions A, B,C, E are defined as in Eq. (15.171) and we set 


G(q) := 


7T 2 — a 2 


127r 2 a 2 p A 


(5.26) 


We notice that, in accordance with the existing literature (see, e.g., HU), the con¬ 
formal part of the renormalized stress-energy VEV coincides with the analogous 
contribution derived for Dirichlet boundary conditions in subsection 15.61 besides, 
comments analogous to those made at the end of the cited subsection also hold in 
this case. Let us only remark that the VEV (0|T /t! ,(q) |0) ren has an additional term 
proportional to the function G(q); this function changes sign for either a < n or 
a > 7T and diverges for p —* 0. 

Concerning the pressure on the boundary, also in this case Eq. (15. 12)1 clearly yields 
a divergent result; on the other hand, using Eq. (15. lip , we obtain 


Pi 


5io 7T 4 — a 4 / 1\ 7T 2 — a 2 

4n 2 p 3 120a 4 \ 6/ a 2 


(5.27) 


As in the previous subsection, we find that the renormalized pressure depends of 
the parameter £. 

Proceeding as explained in subsection 15.51 the renormalized stress-energy VEV 
(I5.25P and pressure (15.27P are easily seen to give for a = tt and a = 7t/ 2, re¬ 
spectively, Eq.s (I4.62p (I4.63p with ot\ = 1, and Eq.s (14.641) (14.651) with ot\ = a 2 = 1. 

5.9 Periodic boundary conditions (the cosmic string). Finally, let us con¬ 
sider the case where the field fulfills periodic boundary conditions on the half-planes 
7r 0 ,7r a , meaning that 


$(t, p, 0, z) = $(t, p, a, z) , d g $(t, p, 0, z) = de^it, p, a, z) 
for t, z G R, p G (0, +oo) . 


(5.28) 
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Ill passing, let us mention that the same framework was also analysed by Dowker 
[[13 ] and by Fulling et al. [23], both employing a point-splitting approach; more 
precisely, in [13] the conformal part of the energy density alone is computed, while 
in [23] the authors only report the graphs of the energy density and pressure for 

C = 1/4. 

Similarly to the cases of the segment and parallel hyperplanes configurations with 
periodic boundary conditions (considered, respectively, in subsection 6.9 of Part I 
and in subsection 13.91 of the present paper), the spatial domain hi for the present 
setting is more properly addressed as a flat Riemannian manifold. The manifold hi 
has a global coordinate system q = (p, 9, z) : 12 — > (0, +oo) x T/ x R, x q(x) 
where the second factor is the one-dimensional torus T* := R/(crZ); the line element 
in these coordinates has the form (15.2ft (0). The corresponding spacetime R x Q 
(with the line element ds 2 = — dt 2 + d£ 2 ) is usually described in terms of a “cosmic 
string” due to the presence of a 1-dimensional topological defect coinciding with the 
axis {p = 0}. 

A complete orthonormal system of (improper) eigenfunctions (F k ) k£ )c of A in L 2 (Q), 
with the related eigenvalues (cuDk&c, is given by 


^(q) := 


co 


J\\n\(up) e 


i\ n 9 gihz 


\r> • 


2nn 


2na a ’ (5.29) 

col := ^ for A; = (n,co,h) e K, = Z x (0, + 00 ) x R 


(similarly to the previous subsections, we are assuming K, to be a measure space 
such that J K dk = J2t=-oo / 0 + °° dco dh). The modified cylinder kernel T can 
then be evaluated starting from its eigenfunction expansion (I2.28P : 


F k (x)F k ( y) = 


1 

27ra 


r e -t^k 

^(t;q,p) = dk - 

Jk. u k 

+°o r+00 

j\ n (6-e') duuJ Xn {u P )J\ n (up') I dh . 

Jo J- 00 Vc o 2 +h 2 


E 


+00 —t\/u' 2 +h 2 

i _ e ifc(*-v) 


(5.30) 


Evaluating the integrals in h and co as in [23] and considering separately the terms 
with positive and negative values of n, we obtain the expression 


+OO 


+OO 


T{ t;q,p) = 


2'Kapp' sinh v 


i+E' 


,—X n v i\ n (9—9‘ 


">+E 


g -i\n(9-9') 


(5.31) 


n =1 


72—1 


15 In other words, Q is a quotient space of the Dowker manifold Poo; this is an infinite-sheeted 
Riemannian surface that can be described in terms of a global coordinate system 

q ; Poo —^ (0, -boo) x R x R , x h-> q(x) = (p(x), 6»(x), 2 (x)) 

(and of the line element (15.21) ). 
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which in turn, summing the geometric series, yields 


r(t;q,p) 


1 sinh(^n) 

27T app' sinh v cosh(^n) — cos(^(0 — 9')) 


(5.32) 


(again, v is defined as in Eq. (15. 15(1 ). Using Eq.s (15.7H5.8jl once more, we find the 
following expression for the renormalized stress-energy VEV: 


A(q) 


/-I 
0 
0 
V 0 


(0\T^{q)\0) ren 


0 

0 

0\ 


/—2 0 

0 

°\ 

1 

0 

0 

/ i\ . . 

0 -1 

0 

0 

0 

-3 P 2 

0 

+ H) G(q) 

0 0 

3 P 2 

0 

0 

0 

V 


\ 0 0 

0 

2 / 


A(q) 


( 27 t ) 4 — a 4 
14407 r 2 a 4 p 4 ’ 


G(q) 


( 27 t ) 2 — a 2 
247r 2 a 2 p 4 


(5.33) 


Let us observe that the above result does not depend explicitly on the angular 
variable 0; this was to be expected due to the homogeneity of the considered con¬ 
figuration with respect to this coordinate. Besides, as for the cases of Dirichlet and 
Neumann boundaries, both the conformal and non-conformal part of the renormal¬ 
ized VEV of the stress-energy tensor diverge near the axis {p = 0}, that is in the 
proximity of the cosmic string. 

In conclusion, we notice that for a = 2ix, in which case the considered configuration 
is equivalent to that of a scalar field on the whole Minkowski spacetime, Eq. (I5.33jl 
gives A(q) = G(q) = 0. So, we have this result with its own interest: when zeta 
regularization is applied to a massless scalar field on the whole (3 + 1)-dimensional 
Minwkowski spacetime, the renormalized stress-energy VEV vanishes identically. 
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